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Abstract.  We  consider  the  problem  of  self-similar  viscous  limits  for  general  systems  of  conservation 
laws.  First,  we  give  conditions  so  that  the  resulting  boundary  value  problem  admits  solutions.  In 
particular  this  covers  the  class  of  symmetric  hyperbolic  systems.  Second,  we  show  that  if  the 
system  is  strictly  hyperbolic  and  the  Riemann  data  are  sufficiently  close  then  the  resulting  family 
of  solutions  is  of  uniformly  bounded  variation  and  oscillation.  Third,  we  construct  solutions  of  the 
Riemann  problem  via  self-similar  viscous  limits  and  study  the  structure  of  the  emerging  solution 
and  the  relation  of  self-similar  viscous  limits  and  shock  profiles.  The  emerging  solution  consists  of 
N  wave  fans  separated  by  constant  states.  Each  wave  fan  is  associated  with  one  of  the  characteristic 
fields  and  consists  of  a  rarefaction,  a  shock,  or  an  alternating  sequence  of  shocks  and  rarefactions 
so  that  each  shock  adjacent  to  a  rarefaction  on  one  side  is  a  contact  discontinuity  on  that  side. 


1.  Introduction 

Consider  the  system  of  conservation  laws  in  one  space  dimension 

(1.1)  dtU  +  dxF{U)  =  0 

where  x  €  IR,  t  >  0,  U(x,t)  takes  values  in  IR^  and  the  flux  function  F  :  IRN  — +  IR^  is 
assumed  smooth.  If  the  matrix  VF(U)  has  real  and  distinct  eigenvalues  then  (1.1)  is  called  strictly 
hyperbolic,  and  its  eigenvalues  (called  characteristic  speeds)  may  be  ordered 

(1.2)  AX(U)  <  A2(U)  <  ...  <  A N(U). 

Let  r\(U) , ... ,  rjv(Z7)  and  l\(U) ,  ... ,  In(U)  be  the  corresponding  right  and  left  eigenvectors.  They 
are  linearly  independent  and  form  a  pair  of  local  bases  in  the  state  space. 

The  Riemann  problem  consists  of  solving  (1.1)  subject  to  data  a  single  jump  discontinuity 

(1'3>  Itl 
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It  describes  the  local  structure  of  BV  solutions  at  points  of  shock  interactions  (DiPerna  [Dp], 
Liu  [Li4])  and  serves  as  a  building  block  for  solving  the  Cauchy  problem  via  the  Glimm  scheme 
(Glimm  [G]).  In  solving  the  Riemann  problem  one  encounters  loss  of  uniqueness  that  has  to  be 
accounted  for  by  imposing  admissibility  restrictions  on  solutions.  For  weak  waves  in  strictly  hy¬ 
perbolic  systems  it  suffices  to  impose  such  restrictions  only  at  shocks.  Lax  [Lai]  in  the  genuinely 
nonlinear  case  and  Liu  [Lii,  L^]  in  the  general  case  provided  comprehensive  shock- admissibility 
criteria  and  obtained  a  unique  solution  of  (1.1),  (1.3)  for  weak  waves.  The  reader  is  referred  to 
Dafermos  [D3]  for  a  thorough  discussion  of  the  issue  of  admissibility.  The  solution  of  the  Riemann 
problem  is  based  on  the  invariance  of  (1.1),  (1.3)  under  dilations  of  the  independent  variables 
(x,t)  — +  ( ax  ,  af),  for  a  >  0.  Because  of  the  expected  uniqueness,  one  seeks  for  solutions  U  =  U(£) 
functions  of  the  single  variable  £  =  j.  The  function  U  is  a  weak  solution  of  the  boundary  value 
problem  (V) 

U'  +  F(U)'  =  0 

(V) 

U(±  00)  =  u± 

subject  to  admissibility  conditions  on  shocks.  The  classical  solution  of  (V)  consists  of  two  steps: 
First  special  solutions  of  rarefaction  waves,  shock  waves  and  contact  discontinuities  are  studied, 
and  are  used  to  construct  the  elementary  wave  curves.  There  is  one  elementary  curve  associated 
with  each  characteristic  field  with  the  parametrization  of  the  curve  serving  as  a  measure  of  the 
strength  of  the  associated  wave.  Second,  it  is  shown  that  the  compound  curves  emanating  from  a 
fixed  U _  give  rise  to  an  invertible  map  that  covers  a  full  neighborhood  of  right  end  states  U+  ( c.f 
[Lai,  LL*]). 

The  objective  of  this  article  is  to  obtain  the  complete  solution  of  the  Riemann  problem  for 
weak  waves  by  an  alternative  approach,  in  the  spirit  of  viscosity  methods.  Namely,  admissible 
solutions  of  (V)  are  constructed  as  e  \  0  limits  of  solutions  to  the  problem  ( Ve ) 

U*  +  F(U)f  =  eU" 

{Ve) 

U (±00)  =  u± , 

with  £  >  0.  The  latter  consists  of  an  elliptic  regularization  of  the  Riemann  operator  in  ( V ).  This 
approach  is  usually  called  self-similar  viscous  limits  and  was  proposed  by  Dafermos  [Di],  who 
motivated  it  by  introducing  an  artificial  ’’viscosity”  regularization  that  preserves  the  invariance 
under  dilations  of  coordinates.  The  study  of  the  Riemann  problem  amounts  to  performing  the 
following  steps: 

(i)  To  construct  solutions  of  the  problem  (7\),  with  e  >  0  fixed. 
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(ii)  To  construct  solutions  of  ( V )  as  e  \  0  limits  of  solutions  of  (Ve). 

(iii)  To  study  the  structure  of  the  emerging  solution. 

Our  interest  in  (Ve)  stems  from  the  connection  with  the  problem  of  viscous  limits.  For  the 
system  of  viscous  conservation  laws 

(1.4)  dtU  +  dxF(U)  =  ed2xU 

subject  to  Riemann  data,  the  invariance  under  dilations  (x,t)  — >  (ax,  at),  a  >  0,  no  longer  holds. 
It  is  a  simple  calculation  to  see  that  the  solution  Ue  of  (1.3  —  1.4)  can  be  expressed  as 

(1.5)  P*(*,0  =  V(|,-j) 
where  F(£,s)  satisfies 

(1.6)  *  (_^  +  F(y)c) 

o 

for  — oo  <  £  <  oo,  —oo  <  s  <  0.  Therefore,  the  viscous  limit  problem  for  Riemann  data  is  a  two 
parameter  problem  and  studying  the  limit  of  U£  as  e  j  0  amounts  to  studying  the  limit  of  V(£,  s) 
as  s  |  0—.  The  problem  (Vc)  arises  when  replacing  the  parabolic  operator  in  (1.6)  by  an  elliptic 
operator;  its  study  is  expected  to  provide  insight  on  the  difficult  problem  of  viscous  limits.  The 
two  regularizations  have  been  compared  for  Burgers’  equation  (Slemrod  [S2]). 

The  notion  of  self- similar  viscous  limits  appears  in  the  articles  [Ka],  [Tui],  [TU2],  [Di].  Tupciev 
[Tui ,  TU2]  used  them  to  formally  motivate  a  shock  admissibility  condition  for  the  Riemann  problem, 
that  amounts  to  the  requirement  that  admissible  shocks  have  associated  shock  profiles.  The  direct 
use  of  self-similar  viscous  limits  is  initiated  by  Dafermos  [Di ,  D2 ]  who  proposed  it  as  an  admissibility 
criterion  and  devised  a  versatile  framework  for  treating  the  analytical  aspects  of  the  problem.  The 
approach  has  been  tried  on  several  examples  of  strictly  hyperbolic  2x2  systems  [Di,  DDp,  KKr, 
STzj,  TZ2],  on  a  system  of  two  equations  that  exhibits  change  of  type  [Si ,  Fa2],  and  on  the  fluid 
dynamic  limit  for  the  Broadwell  model  [STZ2,  Tzi].  It  has  been  established  at  the  level  of  such 
examples  [D2,  Faj,  TZ2]  that  it  dictates  the  same  structure  for  the  Riemann  problem  solution 
as  that  obtained  by  the  usual  shock-admissibility  criteria,  or  by  requesting  that  each  admissible 
shock  has  an  associated  viscous  shock  profile.  In  contrast  to  most  admissibility  criteria,  self¬ 
similar  viscous  limits  penalize  the  whole  wave-fan  simultaneously.  (It  would  be  interesting  to  try 
the  method  in  one  of  the  examples  where  the  requirement  that  each  shock  has  an  associated  profile 
is  not  sufficient  to  guarantee  uniqueness).  Based  on  that  fact,  a  fitting  terminology  would  be  to 
call  admissibility  via  self-similar  viscous  limits  as  the  viscous  wave-fan  admissibility  criterion. 
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Here,  we  pursue  the  method  for  strictly  hyperbolic  systems  of  more  than  two  equations.  We 
address  the  questions  of  existence,  performing  the  e  — *  0  limit,  and  structure  of  the  emerging 
solution.  The  key  step  lies  in  controlling  the  diffusion  induced  wave  interactions  and  obtaining 
uniform  variation  estimates  for  solutions  of  {Ve).  The  article  is  organized  as  follows: 

In  Section  2  we  study  the  question  of  existence  of  solutions  for  (Vc)-  We  show  that  for  any 
system  equipped  with  an  Lp  estimate  the  problem  (Vc)  admits  solutions  for  each  e  >  0.  The 
analysis  applies  to  the  class  of  symmetric  hyperbolic  systems. 

Sections  3  to  7  are  the  core  of  the  article  dealing  with  the  question  of  obtaining  uniform 
variation  estimates  for  families  of  solutions  to  ( Ve ).  Even  for  Riemann  data,  waves  of  different 
families  can  interact  through  diffusion  and  contribute  to  the  total  variation.  Therefore,  one  has  to 
devise  a  scheme  for  measuring  the  variation  of  the  solution  (through  the  individual  waves)  and  to 
calculate  the  effects  of  wave  interactions.  We  refer  to  Section  3,  which  serves  as  an  introduction 
to  this  part,  for  an  outline  of  our  strategy.  The  outcome  is  summarized  in  Theorem  3.1  and  states 
that  if  (1.1)  is  strictly  hyperbolic  and  the  data  U±  are  such  that  |  U+  —  U-  \  is  small,  then  ( Ve )  has 
solutions  that  are  of  uniformly  bounded  and  small  oscillation  and  variation. 

In  Sections  8,  9  and  10  we  solve  the  Riemann  problem  for  strictly  hyperbolic  systems  via  self¬ 
similar  viscous  limits,  under  the  sole  hypothesis  that  \U+  —  U-\  is  small.  The  variation  estimates 
of  Section  7  are  used  in  Section  8  to  establish  the  e  — >  0  limit  and,  more  important,  to  study 
the  structure  of  the  emerging  solution  U  of  (V).  It  turns  out  that  U  consists  of  N  wave  fans 
separated  by  constant  states.  Each  wave  fan  is  associated  with  one  of  the  characteristic  fields  and 
is  either  a  rarefaction,  or  a  shock  satisfying  a  weak  form  of  the  Lax  conditions,  or  a  composite 
wave  consisting  of  an  alternating  sequence  of  shocks  and  rarefactions  so  that  each  shock  adjacent 
to  a  rarefaction  on  one  side  is  a  contact  discontinuity  on  that  side.  In  Section  9  it  is  shown  that 
for  shocks  that  do  not  correspond  to  linearly  degenerate  characteristic  fields  solutions  of  (VE)  have 
the  internal  structure  of  traveling  waves.  Finally,  in  Section  10  we  compare  the  outcome  of  the 
solution  via  self- similar  limits  to  the  classical  solution  of  the  Riemann  problem  [Lai,  Lii,  L^]. 
Both  for  genuinely  nonlinear  systems  as  well  as  for  general  strictly  hyperbolic  systems,  the  same 
structure  results  for  the  Riemann  solution.  The  relation  with  the  Liu  shock  admissibility  criterion 
is  indirect,  and  follows  from  the  fact  that  (a  strict  version  of)  the  Liu  shock-admissibility  criterion 
is  equivalent  to  the  requirement  that  admissible  shocks  have  associated  shock  profiles  (Liu  [Lia], 
Majda  and  Pego  [MP]). 
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2.  Existence  of  Connecting  Trajectories  for  ( Ve ) 

The  objective  of  this  section  is  to  construct  solutions  of  the  problem  (Ve),  for  e  positive  fixed. 
{Ve)  is  a  boundary-value  problem  for  a  system  of  non- autonomous  ordinary  differential  equations. 
First,  it  is  shown  that  L°°  estimates  are  sufficient  to  establish  existence  of  solutions  for  (Ve).  Then  a 
construction  scheme,  originally  proposed  by  Dafermos  [Di],  is  presented  in  Section  2.2.  Existence 
of  connecting  trajectories  then  relies  on  a-priori  estimates,  which  are  established  in  Section  2.3 
under  various  structural  hypotheses  on  (1.1).  Most  notably,  the  analysis  applies  to  the  class  of 
symmetric  hyperbolic  systems. 

2.1.  Preliminaries.  Assume  that  U  is  a  classical  solution  of  ( Vc )  satisfying  the  bound 

(2.1)  sup  |U(£)|<Af. 

*-oo<£<oo 

where  M  is  a  constant  that  may  depend  on  e.  Integrating  the  differential  equation 

(2.2)  e  U"  =  -£!/'  +  F(U)' 
it  is  easily  seen  that  U  satisfies  the  identities 

(2.3)  t/'(O  =  tf'(0)e-£+±e-£  /*  e^VF{U{0)  U'(()  d( , 

£  Jo 

and 

(2.4)  e  U'(()  =  e  U'{ 0)  +  F(U( 0)  -  t  U{ 0  -  F(U(0))  +  f  U(Q  d( . 

Jo 

Using  (2.1),  (2.3)  and  Gronwall’s  inequality,  we  obtain 

(2.5)  \U'(t)\  <  |tf'(0)l  e<2oKI-*2>/2£ , 

where  a  :=  sup  |VF(F)|. 

|V|<Af 

Integrating  (2.3)  over  y/e)  and  performing  a  change  of  variables  in  the  resulting  integrals, 

we  arrive  at  the  identity 

v% 0)  J'  =  -^(P(v^)  -  U(-Ve)) 

(2-6)  +  -F([/(0))  ['  e-fl'it.  -  -  f  F(U(-/i())  d( 

£  J- 1  £  J- 1 

+  ~f  (e«2-^2F{U(y/iO)d<;dt. 

£  J- 1  Jo 
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In  turn,  this  leads  to 


|£f'(0)|  J1  e-«2/2^  <  I  [  y/i \U(VZ)  -  U(-V~e)\ 

(2.7)  +  sup  \F(U(y/i 0)1  (4  +  2  /  f* Ce«*-?W  <%<%) 

-1<£<1  Jo  JO  J 

<-{M+  sup  |F(F)|). 

£  |V|<M 

On  the  other  hand  (2, 1),  (2.4)  and  (2.7)  give 

(2.8)  TOI  <  f  (1  +  ICI)  ■ 

Relations  (2.5),  (2.7)  and  (2.8)  imply  that  any  solution  obeying  the  bound  (2.1)  will  also 
satisfy  the  first  derivative  estimates 

{c  |t|  <  Oa 

f „<«,-<’>/*  |fj;2a  ■  0<£il- 

In  (2.9)  the  constants  C  and  a  depend  only  on  sup  |f/(£)|,  while  the  exponent  becomes 

—oo<£<oo 

negative  for  |£|  >  2a.  In  addition,  (2.2)  yields 

(2.10)  |£i"(e)l<i(«  +  KI)|£!'(OI, 

which  in  conjunction  with  (2.9)  provides  an  estimate  for  the  second  derivatives. 

2.2.  The  Construction  Scheme.  Let  e  €  (0,1]  be  fixed  and  consider  the  two-parameter 
family  of  boundary- value  problems 

-(U'  +  {iF(U)'  =  eU"  -/<£</ 

(2.11) 

U(±l)  = 

with  parameters  /t  6  [0, 1] ,  l  >  1.  The  following  theorem  [Dj,  p.3]  provides  sufficient  conditions 
that  guarantee  the  existence  of  solutions  for  ( Ve )•  We  outline  its  proof  for  the  sake  of  completeness. 

Theorem  2.1  Assume  that  there  is  a  constant  M  depending  at  most  on  U~ ,  U+,  the  function 
F{U)  and  e  ( but  independent  of  p  and  l ),  such  that  any  solution  U(£)  of  (2.11)  satisfies  the  bound 

(2.12)  sup  |[/(£)|<Af. 

-/<£</ 

Then,  there  exists  a  classical  solution  of  (Vc)  denoted  again  by  U(Q  and  defined  on  (— oo  ,oo). 


6 


Proof.  First,  solutions  of  (2.11)  are  constructed  by  means  of  a  continuation  argument.  Given  a 
smooth  function  V  the  solution  W  of  the  boundary- value  problem 


(2.13) 


£  W"(£)  +  £  W'(£)  =  F(V(0)'  -  /  <  £  <  * , 

W(-l)  =  U.  ,  W(+l)  =  U+ 


is  computed  by  the  formula 


(2.14) 


W(  £)  =  U-  + 


J*F(V(C))d( 


a:  rex~^~  F(V(t))  drd( , 


where  the  constant  Uq  €  IR^  is  calculated  by 


U0 


(2.15) 


f_{  e-£d<;  =  (U+-U-)-±  F(V(  0)  d( 

1  fl  ft  T 

+  ~2  I  /  Te  a*  F(V{T))drdC 
£  J-i  Jo 


Set  X  =  C70([— /,/]  ;1Rn)  and 


Sl:={U  eX  :  sup  |tf(£)|  <  M  +  1}  . 

-/<€</ 

X  with  the  sup-norm  is  a  Banach  space  and  fl  is  a  bounded,  open  subset  of  X.  Consider  the 
map  T  :  -»  X  carrying  V  €  fl  to  W  =  T(V)  defined  by  the  relations  (2.14)  and  (2.15).  T  is 

compact  and  continuous,  and  classical  solutions  of  (2.11)  are  identified  with  fixed  points  of  pT. 
The  map  I  -  pT  :  fi  x  [0, 1]  -*■  X  satisfies  the  hypotheses  of  the  Schaeffer  fixed-point  theorem  (e.g. 
Rabinowitz  [R,  Ch  V]).  Hence,  for  each  p  6  (0, 1]  there  is  at  least  one  solution  of  the  equation 
U  —  pT(U)  =  0  in  the  set  fl. 

Let  now  £/£•;/)  denote  a  solution  of  (2.11)  for  p  =  1.  In  the  last  step,  solutions  of  (Vc)  are 
constructed  as  /  — ►  oo  limits  of  U(- ;  l).  Proceeding  as  in  the  derivation  of  (2.9)  and  (2.10),  it 
follows  that  such  solutions  satisfy  the  bounds  (2.12)  and 

|f/'(£ ;/)|  <  -  e(2«|£|-£2)/2C 

(2-16)  £ 

\u"(£ ;  01  <  ^7  (1  +  |£|)  e(2a|$|-€  )/2e 

with  C  and  a  depending  on  M  but  not  on  l.  Extend  U(-;l)  outside  [—/,/]  by  setting  l/(£;/)  =  li¬ 
ter  £  <  — /  and  [/(£;/)  =  U+  for  £  >  /.  The  Ascoli-Arzela  theorem,  together  with  a  diagonal 
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argument,  implies  the  existence  of  a  sequence  {/„},  ln  00 ■>  and  a  function  U  6  C1((— 00, 00) ;  IR^) 
such  that  >  U  and  U'(- ; ln)  — >  U'  uniformly  on  compact  subsets  of  IR.  Because  of  (2.16) 

the  convergence  is  uniform  on  IR  and  f/(±oo)  =  U±.  Passing  to  the  limit  ln  —>  00  shows  that  U  is 
a  classical  solution  of  (Ve).  | 

2.3  The  a-priori  estimates.  The  scope  of  this  section  is  to  provide  the  sup-norm  estimates 
that  authorize  application  of  Theorem  2.1.  In  the  sequel,  U (£)  stands  for  a  solution  of  the  family 
of  boundary- value  problems  (2.11)  defined  on  [—/,/]  and  depending  implicitly  on  g, ,  l  and  £.  In  the 
process  of  estimating  U (£)  we  pursue  ideas  that  were  developed  by  Dafermos  and  DiPerna  [DDp] 
in  the  context  of  2  x  2  systems  and  use  the  concept  of  entropy-entropy  flux  pairs  (Lax  [La2]). 

A  scalar-valued  function  q(U)  is  called  an  entropy  for  (1.1),  with  corresponding  entropy  flux 
q(U),  if  every  smooth  solution  satisfies  the  additional  conservation  law 

(2.17)  dtv(U)  +  dxq(U)  =  0 

Such  pairs  ( q(U),q(U ))  are  generated  by  solving  the  system  of  (linear)  differential  equations 

(2.18)  Vq(U)  =  Vq(U)  VF(U) . 

Trivial  examples  of  solutions  are  provided  by  (c  •  U  ,  c  ■  F(U)),  with  c  any  constant  vector  in  IRN. 
Since  (2.18)  is  overdetermined  for  N  >  3,  for  systems  of  three  or  more  equations  the  existence 
of  (nontrivial)  entropies  is  the  exception  rather  than  the  rule.  Nevertheless,  specific  systems  that 
arise  in  applications  are  often  naturally  endowed  with  some  entropy-entropy  flux  pairs.  Also,  the 
class  of  symmetric  hyperbolic  systems,  that  is  systems  for  which  VF(U)  is  a  symmetric  matrix, 
admits  the  pair 

(2.19)  V(U)  =  \  \U\2  q(U )  =  U  •  F(U)  -  g(U) , 

where  g  is  a  potential  for  F  satisfying  F(U)  =  Vg({7). 

Let  ( r](U),q(U ))  be  an  entropy-entropy  flux  pair  for  (1.1).  Using  (2.18)  we  deduce  that 
solutions  of  (2.11)  satisfy  the  identity 

(2.20)  -tq'  +  fiq'  =  er]"  -  eU'  -(V2t])U' 

where  77  =  v{U(0)  >  9  =  q(U(0)-  exploiting  (2.20),  it  is  helpful  to  use  entropy  functions 
77(1/)  that  are  convex  (or  linear).  The  following  lemma  indicates  how  to  bound  the  total  entropy 
production.  Given  a  constant  entropy  level  77,  consider  the  level  set 

(2.21)  Cf,  =  {U  £MN  :  q(U)=fj}. 
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If  Cf)  is  nonempty,  let 

(2.22)  Qn=  sup  \q(Ux)-q(U2)\ 

Vi  ,t/2€C„ 

be  the  oscillation  of  q(U)  on  the  level  set  Cjj. 

Lemma  2.2.  Assume  that  T)(U)  is  a  convex  entropy  with  corresponding  entropy  flux  q(U).  If  fj  is 
any  constant  such  that 


(2.23)  r]  >  max  {v(pU-), r)(pU+)} 

0<M<1 

then  for  any  (a,/?)  C  (-1,1) 

(2.24)  [P  (V(U(Z))-v)d{<K, 

J  a 

where  K  —  if  rj(U(£))  >  fj  for  some  f  €  (a,/?),  and  K  =  0  otherwise. 

Proof.  The  proof  is  based  on  the  following  observation.  Let  fj  be  a  fixed  entropy  level  and  suppose 
that  a ,  b  are  two  points  in  (—1,1)  with  the  properties  a  <  b  and 

(2.25)  v(U(a))  =  T](U(b))  =  fj  with  (rj  o  U)'(a)  >  0 ,  (i)  o  U)'(b)  <  0 . 


Integrating  (2.20)  over  [a,  6],  we  obtain 


(2.26) 


I*  iv(U(0)  -f))d(  +  s  f  U'( 0  •  V2v(U(0)  U'(0  d£ 

J  a  J  a 

<  -p[q(U(b))-q(U(a))]  <  Qn, 


which,  upon  using  the  convexity  of  f](U),  yields 


(2.27) 


[b  {n(u(0)-v)d(i<  Qf). 

J  a 


If  y(U(0)  <  fj  for  —l  <  f  <  l,  then  (2.24)  is  trivially  true  with  K  =  0.  So  suppose  that  the 
set  {£  €  (—1,1)  '■  v(U(£))  >  fj}  is  nonempty.  It  is  also  open  and  thus  admits  a  decomposition  into 
a  countable  union  of  disjoint  subintervals 


(2-28)  {£  e  (-1,  l) :  V(U(0)  >  V}  =  IJ  (°*> M > 

*€/ 

where  k  ranges  over  an  index  set  I  (either  a  finite  set  or  the  integers).  For  fj  restricted  by  (2.23) 
the  points  a*  and  6*  lie  in  (-1,1).  Also,  since  tj(U( £))  >  fj  for  a*  <  f  <  bk  with  k  €  I,  relations 
(2.25)  are  satisfied  at  the  endpoints  a k  ,  bk. 
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Given  any  (a,/?)  C  (—1,1),  choose  a,  b  as  follows:  If  r)(U(a ))  >  fj  set  a  =  sup{a*;  <  a},  while 
if  Tj(U(a))  <  fj  set  a  =  inf{afc  >  a};  if  T}(U(fi))  >  fj  set  b  =  inf{f>fc  >  /?},  while  if  r](U(fi))  <  fj  set 
b  =  sup{6jt  <  /?}.  If  q(U(£))  >  fj  at  some  £  €  (ot,fi),  a  and  b  are  well  defined,  a  <  b,  relations 
(2.25)  are  satisfied  at  a,  b  and 

(2.29)  f  (V(U( {))  -  <  f  (mo>  -?)«<«»• 

*/a 

Otherwise  (2.24)  holds  with  A  =  0.  | 

In  general  the  quantity  depends  on  the  form  of  the  level  set  as  well  as  the  function 
q({7)  and  may  be  infinite.  If  it  happens  that  is  a  compact  set,  then  is  finite  and  (2.24) 
provides  an  integral  estimate  independent  of  p  ,1  and  e.  An  entropy  is  called  normal  if  ij(U)  — *  oo 
as  |£/|  — y  oo.  If  the  system  (1.1)  is  endowed  with  a  convex  normal  entropy,  then  nonempty  level 
sets  Cfj  are  compact,  and  that  leads  to  integral  estimates  of  the  type  (2.24).  For  a  symmetric 
hyperbolic  system  ij(U)  =  \\U\2  is  an  example  of  a  convex  normal  entropy. 

Next,  we  present  two  approaches  for  obtaining  the  sup-norm  estimates  (2.12).  The  first  exploits 
the  entropy  identity  (2.20),  and  requires  the  existence  of  a  strictly  convex,  normal  entropy  function 
t i(U),  defined  (only)  on  the  exterior  of  some  open  ball  in  the  state  space. 

Proposition  2.3.  Assume  that  (1.1)  admits  a  strictly  convex,  normal  entropy  tj(U )  defined  on 
the  exterior  of  a  ball  and  satisfying  the  growth  restriction:  There  are  q  >  0  and  positive  constants 
C  and  ro  such  that 

(H)  |V7?(t/)|2  <Cv(U)rj(Uf-q  for  |J7|>r0, 

where  v(U)  is  the  smallest  eigenvalue  of  the  Hessian  V2t/(1/).  Then  solutions  of  (Ve)  exist  for 
every  e  >  0. 

Proof.  Let  ri(U)  be  a  strictly  convex,  normal  entropy  defined  for  {17  €  IR^  :  |17|  >  ro}  and 
satisfying  (H)  for  some  q  >  0.  Without  loss  of  generality  we  may  assume  that  rj(U)  is  positive. 
Let  U(£)  be  a  solution  of  (2.11)  on  (—1,1).  For  those  £  that  |tl(£)l  >  ro  equation  (2.20)  is  satisfied. 

Let  r  >  max{|17+| ,  |(/_| ,  ro}  and  fjr  =  max|(/|=r  q(U)  be  fixed  and  choose  two  entropy  levels 
fji  >  f)i  >  f)T  >  0.  Consider  the  set 

(2.30)  A  =  {i  €  (-/,  /)  :  rj(U( 0)  >  th  ,  \U(0\  >  r  }  . 

Since  tj(U)  — ►  oo  as  \U\  — ►  oo,  if  the  set  A  is  empty  then  sup_/<^<;  \U(£)\  <  M,  for  some  M 
depending  on  fji  and  r,  and  thus  (2.12)  holds  in  this  case.  So,  assume  that  A  is  nonempty.  It  is 
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also  open  and  thus  admits  the  decomposition  A  =  Ujfce/(a*>  &*)  into  a  countable  (or  finite)  union 
of  disjoint  intervals.  In  addition  the  choice  772  >  fjr  implies  that,  for  any  k  G  /, 

T){U  (0)  >  j?2  for  ak<Z<bk, 

(2.31) 

v(U(ak))  =  rj(U(bk))  =  1)2  and  (77  o  U)\ak)  >  0 ,  (v  °  U)'(bk)  <  0 . 

Henceforth  we  focus  on  a  fixed  interval  ( ak,bk ).  Let  rk  be  a  point  where  p(U(£))  assumes 
its  maximum  in  the  closed  interval  [ak,bk].  Using  Schwarz’s  inequality,  the  strict  convexity  of  77, 
hypothesis  ( H )  and  relations  (2.31),  (2.25)  and  (2.26)  we  obtain 

(2/?)  [,(t?(rt))«  -  =  r  VijW0)P'«K 

Jak 

(2.32)  <  ( £  S  {£  U'(0  -  V2^<U(QjU'(C)d(\ 1 

<  [C £  ■ 

For  those  U  G  IRN  that  p{U )  >  >  fji  >  0,  it  is 

(2.33) 

Then  (2.32)  yields  the  estimate 

(2.34)  <  (th)i  +  (g/2)  ( 0«.  )’/a  ( ( r>(tf(C))  -  m  )  dC ) * . 

Set  a  =  inf{£  G  (-/,afc)  :  if(ff(0)  >  %on(^,afc)},  6  =  sup{£  G  (&*,/) :  J?(U(C))  >  7/1  on  (&fc,0>- 
Since  |f/(±/)|  <  r  and  771  >  fjr ,  a,  6  are  well  defined  and  satisfy  —l<a<ak<bk<b<l.  In 
addition  (2.25)  holds  and,  as  in  the  proof  of  Lemma  2.2, 

(2.35)  f*  (77(17(0)  -  m)  d{  <  f  (v(U( 0)  -  m)  dt  <  Qni . 

J  ak  J  a 

As  a  consequence,  the  right  hand  side  of  (2.34)  is  bounded  independently  of  k,  and  (2.12)  holds  in 
the  case  that  A  is  nonempty  too.  The  conclusion  now  follows  from  Theorem  2.1.  | 

Regarding  the  growth  assumption  (H)  the  following  remarks  are  in  order.  If  the  strictly 
convex,  normal  entropy  function  is  of  the  form  r)(U)  =  (1/p)  \U\P,  with  p  >  1,  one  easily  calculates 

(2.36)  Vt7(U)  =  \U\P~2  U ,  V2t7(U)=  \U\P~2 1  +  (p  -  2)\U\p-4  U  ®  U . 

The  Hessian  of  77  is  a  positive  definite  matrix  having  eigenvalues:  (p—  l)|t/|p-2  with  corresponding 
eigenvector  U ,  and  \U\P~2  of  multiplicity  N  —  1  with  corresponding  eigenvectors  U 1  any  vector 
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orthogonal  to  U.  Hypothesis  (II)  is  then  satisfied  with  q  =  2.  On  the  other  hand,  if  tj(U)  grows 
like  a  power  up  to  first  order  derivatives,  i.e., 

(2.37)  -c\U\p  <  i 7(17)  <  c\U\* ,  |V»?(C/)|  <  c\U\*~x  , 

for  some  positive  constant  c,  then  ( H )  becomes  a  restriction  on  the  decay  bf  the  minimum  eigenvalue 
for  U  large  and  is  satisfied  provided  that  v(U)  >  \U\~3  for  some  s  <  p  +  2. 

As  a  consequence  of  the  above  remarks  in  conjunction  with  Proposition  2.3,  we  have. 

Theorem  2.4.  If  (1.1)  is  a  symmetric  hyperbolic  system,  then  solutions  of  (Ve)  exist  for  every 
£  >  0. 

In  the  interest  of  developing  technique,  we  present  an  alternative  way  for  establishing  (2.12) 
for  symmetric  hyperbolic  systems.  The  actual  result  is  weaker  than  Theorem  2.4,  as  it  requires  a 
growth  assumption  on  the  flux  F(U),  but  the  approach  may  be  useful  for  other  problems. 

Proposition  2.5.  Suppose  that  (1.1)  is  a  symmetric  hyperbolic  system,  such  that  the  flux  function 
satisfies  the  growth  assumption 

(2.38)  \F(U)\  <  C  (1  +  \U\)P 

for  some  positive  constants  C  and  p  <  3.  Then  solutions  of  (Ve)  exist  for  every  £  >  0. 

Proof.  Symmetric  hyperbolic  systems  are  endowed  with  the  entropy  -  entropy  flux  pair  (2.19),  for 
which  (2.20)  takes  the  form 

(2.39)  -f  (\U\2)'  +  2/x  (U  •  F(U)  -  g(U))'  =  £  (\U\2)"  -  2e \U'\2 

The  function  g  is  a  potential  for  F  satisfying  F(U)  =  Vg(U).  It  can  be  defined  by  the  formula 

(2.40)  g(U)  =  £  jtg(t  U)  dt  =  jf '  F(t  U)  ■  V  dt , 

where  g  has  been  normalized  by  setting  <jr(0)  =  0.  Assumption  (2.38)  induces  a  growth  restriction 
on  g  as  follows: 

(2.41)  \g(U)\  =  |  f1  F(tU)  -Udt\<C(l+  \U\)P+1 . 

Jo 

Set  r  =  max{|{7_|, |f/+|}  and  consider  any  point  f  €  (—/,/)  such  that  |!7(f)|  >  r  and 

(d\U\2/d£)(£)  >  0.  Define  £'  =  inf{C  €  (£,/]  :  |£f(C)|  <  1^(01)  and  observe  that  is  well 
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defined  with  £  <  <  /.  Moreover,  |tf(f')l  =  \U(t)\,  (d\U\2 /dZ)(t')  <  0  and  \U(()\  >  |?7(^)|  for 

£  <  C  <  f  Integrating  (2.39)  over  [£,  £'],  we  obtain 

-/  <(|tf)'«MC  +  2£  /  l^'tOI^C 

+  2„  [  If«')  ■  F(tf({’))  -  sWO)  -  tf(f)  ■  W(0)  +  »(»({))] 

=  -  e(rf|i/l2/^)«) 

-J  C(\u\2)'(Od(  =  (|i'(OI!-|£'({)l2)rf(>  0, 

(2.42)  together  with  (2.38)  and  (2.41)  yield 

(2.44)  e®(0  <  8C(1  +  \Ump+1- 

Note  that  the  bound  (2.44)  holds  for  any  £  €  (— 1,1)  such  that  |I7(f)|  >  r. 

To  conclude  the  proof  fix  two  levels  r\  and  r2,  with  r2  >  rj  >  r,  and  consider  the  set 
B  =  {£  G  (-/,/)  :  |tf(0|  >  r2 }.  If  5  is  empty  then  (2.12)  holds  and  Theorem  2.1  implies  the 
desired  result.  If  B  is  nonempty,  then  it  can  be  decomposed  into  an  at  most  countable  union  of 
disjoint  subintervals  (a*, 6*)  such  that  |tl(a*,)|  =  |t/(£>fc)|  =  r2  and  |I7(f)|  >  r2  for  a*  <  f  <  bk. 
In  each  of  the  intervals  [<ifc,fefc]  the  differential  inequality  (2.44)  is  satisfied.  Next,  fix  k  and  let 
Tk  €  [ofe,hfc]  be  a  point  where  |f/(rfc)|  =  maxot<{<6t  |tf(f)|.  Lemma  2.2,  applied  for  the  entropy 
T](U)  =  \U\2  and  the  level  f)  =  r2,  implies 

(2.45)  fV«)l2-<-M<<3,;<oo. 

Jak 

Since  the  ratio  |i/|1-p(l  +  \U\)1+P/(\U\2  -  r2)  remains  bounded  for  |t/|  >  r2,  using  (2.44)  and 

(2.45)  we  deduce 

(2.46)  e  £  MOl1-”  if  <  C'£  (TOP  -  r2)  d{  <  C'Qrj . 

In  turn,  performing  the  integration  in  (2.46)  yields 

(2-47)  \U(Tk)\3~p  <  (r2 )3_p  +  , 

for  p  <  3,  and 

(2.48)  \U(rk)\<r■leC'Q•\n, , 

for  p  =  3.  In  either  case  (2.12)  holds  and  proof  is  complete.  | 


(2.42) 

Since 

(2.43) 
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3.  Solution  decomposition  -  The  main  result 

The  aim  of  this  article  is  to  construct  solutions  of  the  Riemann  problem  (V)  as  £  —*■  0  limits 
of  solutions  of  ( Ve ).  The  central  difficulty  lies  in  obtaining  f-independent  variation  estimates  for 
families  of  solutions  of  the  problems  (' Pe ).  The  reason  is  that,  even  for  Riemann  data,  there 
are  wave  interactions  induced  by  the  coupling  through  the  self-similar  viscosity  that  need  to  be 
accounted  for.  The  derivation  of  the  variation  estimates  follows  from  a  lengthy  analysis,  carried 
out  in  Sections  3-7.  The  present  section  serves  as  an  introduction,  where  we  outline  the  general 
strategy,  introduce  the  main  hypotheses,  and  present  certain  interesting  geometric  properties. 

Our  approach  is  motivated  by  a  detailed  study  of  the  following  question:  Suppose  we  are  given 
a  family  of  solutions  to  ( VE )  of  uniformly  small  oscillation 

(Co)  sup  \Ue(t)  -U.\<fi. 

—  oo<£<oo 

Such  a  family  would  also  satisfy  uniform  L°°  bounds 
(Cb)  sup  \Ue(0\<M, 

-oo<«oo 

where  the  constants  M  and  /(  are  independent  of  e  and  \i  is  also  small.  Examine  under  what 
structural  hypotheses  on  (1.1)  the  family  {U£}c>o  is  of  uniformly  bounded  variation 


(■ S )  TV^^U'KC. 

It  is  instructive  to  give  a  proof  of  ( S )  for  the  single  conservation  law,  which  contains  some 
ingredients  of  the  approach  followed  for  systems.  Let  {u£}£>o  be  a  family  of  scalar- valued  function 
satisfying 

eu"  =  -£< +  /(«*)' 

(3.1) 

ue(±oo)  =  u±  . 

and  the  uniform  bounds  (Cb)  (which  are  easily  justifiable  in  this  case).  Let  A(u)  =  f'(u)  be  the 
characteristic  speed  of  the  associated  hyperbolic  equation.  It  is  easy  to  see  that  solutions  of  (3.1) 
satisfy  the  representation  formula 

(3-2)  n'M)  =  («+  -  »-)  f 


where 


From  the  form  of  (3.2),  it  follows  that  {u^}  are  uniformly  bounded  in  L1  and  thus  {ue}  is  of 
uniformly  bounded  variation. 

Returning  to  the  general  case,  we  note  that  the  system  (1.1)  is  assumed  strictly  hyperbolic, 
but  no  other  structural  assumptions  are  imposed.  The  eigenvalues  of  VF(U)  are  denoted  by 

(3.4)  Ai (U)  <  A 2(U)  <  ...  <  XN(U) 

and  are  ordered.  The  corresponding  right  eigenvectors  ri(U) , ... ,  r^(U)  and  left  eigenvectors 
l\(U) ,  ... ,  lN(U)  are  linearly  independent  and  satisfy  the  relations 


(3.5) 

VF(U)ri(U)  =  Xi(U)ri(U), 

(3.6) 

li(U)-VF(U)  =  Xi(U)li(U), 

(3.7) 

wwid-{J„  it] 

{r<}  and  {/,•}  form  a  pair  of  local  bases  in  the  state  space  IRW.  By  normalizing  one  of  these  bases 
we  can  attain 

(3.8)  U(U)-rj(U)  =  Sij. 

The  family  {U£}c>o  consists  of  solutions  to  the  boundary  value  problem  (Ve)  that  connect 
two  fixed  end  states  U-  and  U+.  Conditions  that  guarantee  existence  of  solutions  for  (Vc)  are 
given  in  Section  2;  nevertheless,  the  forthcoming  analysis  is  independent  of  such  considerations, 
and  eventually  it  will  also  suggest  a  construction  scheme.  We  assume  the  members  of  {Ue}e>o 
satisfy  the  hypothesis  (C0)  of  uniformly  in  e  small  oscillation  and  (a-fortiori)  the  uniform  bound 
(Cb).  That  restricts  the  data  U±  to  satisfy 

( Hd )  \U+  —  U~\  <  t 

with  r  sufficiently  small.  Also,  each  wave  speed  is  bounded 

(3-9)  Xk- <  Xk{UE(0)  <  **+ 

by  constants  Ajt_,  A^+  independent  of  e.  By  choosing  fi  sufficiently  small,  we  guarantee  that  the 
wave  speeds  are  totally  separated  along  the  family  {f7e}e>o  i  that  is 

Ai-  <  X^U'iO)  <  <  A2_  <  x2{u€(0)  <  A2+  <  ... 

(3.10) 

<  <  Aw— 1  (f^e(O)  ^  ^(N-X)+  <  A N-  <  X n{Uc(0)  ^  ^ N+  ■ 
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The  bound  ( C &)  implies  the  derivatives  of  U£  satisfy  the  estimates  (2.9)  and  (2.10)  with  the  con¬ 
stants  C  and  a  independent  of  s.  In  the  sequel  we  use  the  following  conventions  on  notation  : 
The  ^-dependence  will  be  suppressed  from  functions,  except  at  places  where  emphasis  is  needed. 
By  contrast,  any  ^-dependence  of  constants  will  be  explicitly  stated  by  either  recording  the  precise 
dependence  or  by  using  £  as  a  subscript. 

Consider  the  decomposition  of  U'e  in  the  basis  of  right  eigenvectors  evaluated  at  the  local  value 
of  the  solution  U£ 

N 

(3.H)  ^(0  =  Z>(0r*(tf.(0)- 

k= 1 

The  amplitudes  a k  can  be  recovered  by  using  (3.8) 

(3.12)  ofc(0  =  lk(U£(0)-U'c(0- 


Also,  integrating  (3.11)  over  (—00,00),  we  have 

N  ,00 

(3.13)  U+-V-  =  Y,  /  *k(Qrk(U'( 0)d<- 

k= 1  J~°° 

To  compute  the  equations  that  a*  satisfy,  take  the  inner  product  of  (2.2)  with  lk(U£)  to  obtain 


(3.14) 
and  hence 


-tak  +  A k(U£(0)ak  =  elk(U£(0)  •  U? 

=  ea'k-eVlk(Ue(0)U'e-U'e , 


N  N 

(3.15)  £  a'k  +  [£- A*(J7«(0)]  =  £  E  E  [V/*(CT«(0)  rm(U£(0)  •  rn(U.(0)]  am  an 

m= 1  n— 1 


If  we  introduce  the  notation 


(3.16) 


At  =  A  k(U£(0) 


(3.17)  (3k<mn  =  0k,mn(U£( 0)  =  V/fc(£fe(0)  rm(Ue(0) '  rn(Ue( 0) 

then  a,k  satisfy  the  coupled  system  of  ordinary  differential  equations  with  variable  coefficients 

N  N 

(3.18)  e  afk  +  (£  —  A*)  a*  =  £  ,mn  • 

m=l  n— 1 
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At  this  point  several  remarks  are  in  order.  First,  the  decomposition  (3.11)  is  partly  motivated 
by  the  classical  solution  of  the  Riemann  problem  (Lax  [Lai],  Liu  [L^]).  It  is  expected  to  capture  the 
behavior  near  rarefactions,  but  it  is  not  a-priori  clear  that  it  should  work  well  near  shocks.  Good 
overall  performance  would  indicate  that  (3.11)  captures  the  nature  of  diffusion  induced  averaging 
at  a  shock.  The  coefficients  fik<mn  depend  on  the  solution  and  on  geometric  characteristics  of  the 
surface  S  =  F(U).  The  quadratic  terms  in  (3.18)  represent  the  effect  induced  on  the  Ar-family  by 
interactions  of  waves  of  all  the  families,  and  Pk>mn  measure  the  weights  of  such  contributions.  By 
virtue  of  (Cf,),  Pk,mn  are  uniformly  bounded 

(3-19)  \/3k>mn\  <B. 

Let  gk  be  the  antiderivative  of 


(3.20) 


g'k  =  t-\k  =  t;-\k(Uc(0) 


defined  within  an  arbitrary  constant  of  integration  by 
(3.21) 

In  view  of  (3.9),  it  is 


9k  =  f  s-  Xk(U£(s))ds 
J  a 


(3.22) 


s-  Xk+  <  s-  A k(U£(s))  <  s  -  Afc_  , 


which  in  turn  implies  g'k  >  0  for  £  >  Afc+  ,  g'k  <  0  for  £  <  A*.,  and  gk  looks  like  a  potential- 
well  function  (see  Figure  1).  Let  pkc  be  a  point  where  gk  attains  its  global  minimum,  gk(pkc )  = 
min  </a:(€)-  Then  A*_  <  pkc  <  Xk+  while  the  value  of  gk(pk£)  depends  on  the  choice  of  the  arbitrary 
constant  in  (3.21).  By  setting 

(3-23)  </*(£):=  s-\k(U£(s))ds 

J  Pkc 

we  attain  gk(£)  >  gk(pkc)  =  0  for  £  6  IR.  Furthermore,  Xk(Uc(pke))  =  pke  and  gk(£)  =  0(|£|2)  as 

Ifl  00 • 

Consider  the  linearization  of  the  system  (3.18).  It  consists  of  the  decoupled  system  of  equations 
(3.24)  e  <p'k  +  (f  -  Afc)  (pk  =  0 , 


whose  solutions  are  constant  multiples  of 


(3.25) 


<Pk  = 


e~\9k 


f- 


°°  e~igt  d( 


-7  T  s-*k(Ut(s))  ds 


foo  -if  s-\k(Ut(s))ds 

/  e  Jpkt 

J  —  oo 
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Due  to  their  form  {<pke}  are  strictly  positive  functions  that  are  uniformly  (in  £)  bounded  in  L1 . 


A  comparison  of  (3.24)  with  (3.1)  and  (3.25)  with  the  representation  formula  (3.2)  shows  that, 
for  the  case  of  the  single  equation,  it  is  precisely  the  above  step  that  provides  the  variation  bounds. 
Due  to  the  quadratic  terms  in  (3.18)  though,  this  is  insufficient  for  systems  of  conservation  laws. 
There  are  two  problems  that  we  need  to  account  for,  in  the  case  of  systems.  First  to  understand 
the  effect  of  the  quadratic  terms.  Second,  differential  systems  like  (3.18)  are  best  handled  with 
pointwise  conditions.  On  the  other  hand  the  only  existing  information  (3.13),  relating  the  data  U± 
with  the  amplitudes  a k,  is  of  integral  type.  It  is  thus  necessary  to  devise  a  scheme  that  connects 
pointwise  with  integral  information. 

We  proceed  by  introducing  a  decomposition  of  a*  in  the  form 

(3.26)  ak  =  rk  <pk  +  0k  , 
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where  <pk  is  given  by  (3.25)  and  0k  satisfies  the  system  of  differential  equations 

N  N 

(3.27)  £0>k  +  (t-xk)0k  =  e  EE*  ,mn  {Tm  tym  +  ^m)  (^n  Vn  +  $n)  • 

m=l  n=l 

Clearly  the  sum  Tk  <pk  +  6k  is  a  solution  of  (3.18).  The  idea  is  to  seek  an  asymptotic  expansion 
of  the  wave  amplitude  a*  in  a  parameter  r  =  (tj,  ... ,  tjv),  where  rk  is  thought  as  a  measure  of 
the  strength  of  the  &-th  wave,  and  to  construct  an  expansion  uniform  in  e  in  the  i1-norm.  In 
this  expansion  Tk<pk  is  the  leading  term  and  6k  is  the  error,  which  should  be  of  order  0(|r|2)  as 
|r|  =  |ri|  +  ...  +  \tn\  —*  0.  Clearly,  such  an  expansion  depends  on  the  provided  data,  and  the  key 
question  becomes  under  what  conditions  to  solve  (3.27). 

Next,  we  outline  the  strategy  we  follow  and  the  attained  results  concerning  those  problems  : 
Fix  ci,  C2,  ...  ,cn  to  be  the  respective  middle  points  of  the  intervals  [Ai_,Ai+],  [A2_,A2+],  ...  , 
[A;v_,  A^r+].  Given  a  constant  vector  r  =  (tj,  T2,  ... ,  r>/)  E  IR‘v,  we  consider  (3.27)  subject  to  the 
conditions 

(3.28)  6k{ck)  =  0 , 

and  for  |r|  sufficiently  small  construct  a  solution  6k(£\  T)  that  satisfies  the  estimate 

N 

(3.29)  |0*(-;r)|<C|r|3 

m=l 

This  construction  is  performed  in  Section  5.  It  is  based  on  detailed  estimates,  that  are  presented 
in  Section  4,  on  the  functions  (pk  and  on  integrals  involving  (pmipn  and  capturing  wave  interac¬ 
tions.  The  method  is  to  apply  the  uniform  contraction  principle  to  a  weighted  space  of  continuous 
functions.  The  selection  of  the  weight  is  motivated  by  the  analysis  of  Section  4.  The  analysis  of 
Section  5  validates  the  asymptotic  expansion 

(3.30)  ak(- ;  r)  =  rk  <pk(-)  +  0k(- ;  r) 

for  the  amplitude  a*  in  the  parameter  r.  Note  that  a*  satisfy  the  pointwise  information 

(3.31)  ak(ck ;  r)  =  r*  <pk(ck) 
and  solve  (3.18)  but  not  necessarily  (3.13). 

The  objective  of  Section  6  is  then  to  show  there  exists  a  choice  of  r  =  (tj,  ... ,  r n)  such  that 
(3.13)  is  fulfilled.  To  this  end,  we  consider  the  map  S  :  IR^  — ►  JHN  that  connects  the  wave 
strengths  to  the  boundary  data  by  taking  r  to 

N  fOO 

(3-32)  s(r)  =  /  K  V’fc(C)  +  T)]  rk(Ue(0) • 

A=1  J~°° 
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We  show  in  Section  6  that  S  is  locally  invertible  in  a  neighborhood  of  r  =  0,  and  that  the  inverse 
map  S~l  is  uniformly  bounded  independent  of  e. 

In  Sections  4  to  6,  we  identify  the  precise  hypotheses  (supplementary  to  (3.10)),  on  the  behavior 
of  the  wave  speeds  A*  and  the  right  and  left  eigenvectors  rk  and  lk  along  solutions  U£,  that  are 
necessary  to  carry  out  the  intermediate  steps.  All  these  hypotheses  are  fulfilled  if  the  oscillation 
of  the  family  {£4}£> o  is  restricted,  uniformly  in  e.  It  is  convenient  to  phrase  the  analysis  by 
using  a  general  function  V  of  restricted  oscillation,  sup^eIR  |F(£)  —  U-\  <  p,  in  the  place  of  a 
member  of  {[/£}£>o  Apart  from  splitting  naturally  the  various  parts  of  the  analysis,  this  has 
another  advantage.  The  considerations  of  Sections  4  to  6  motivate  a  construction  scheme  that 
enables  us,  given  Riemann  data  U±  with  \U+  —  U-  \  small,  to  use  the  Schauder  fixed  point  theorem 
and  construct  solutions  Ue  of  (V£)  that  are  of  uniformly  small  oscillation  as  well  as  of  uniformly 
small  variation.  One  interesting  feature  of  the  scheme  is  that  it  is  based  on  the  quadratic  equation 
(3.18)  rather  than  on  a  linearized  equation.  This  final  part  of  the  analysis  is  carried  out  in  Section 
7.  It  justifies  in  particular  Hypothesis  ( Ca )  and  leads  to  the  following  theorem. 

Theorem  3.1  Assume  that  (1.1)  strictly  hyperbolic  and  let  U~  be  fixed.  There  exists  r  sufficiently 
small  such  that  for  e  >  0  and  \U+  —  U-  |  <  r  the  problem  (V£)  has  a  solution  Ue  with  the  properties: 

(i)  The  family  {Ue}e>o  satisfies  (C0)  with  some  fi  independent  of  e. 

(ii)  The  solutions  U£  satisfy  the  representation  formula 

N 

(3.33)  U'£  =  £  [  Tki£<pk  +  <?*(•;  re)}  rk(U£) , 

k= 1 

where  <pk  is  given  by  (3.25),  0*(-;r)  satisfies  (3.29),  and  re  solves  S(re)  =  U+. 

(iii)  The  family  {Uf}c>o  is  uniformly  bounded  in  Z1(IR)  and  {U£}c> o  is  of  uniformly  bounded  (and 
small)  variation. 

We  list  below  certain  geometric  properties  of  the  surface  S  =  F(U)  relating  to  the  coefficients 
fik.mn  =  V/i  rTO  •  rn.  First,  {rk}  and  {lk}  form  bases  of  the  tangent  and  cotangent  spaces  of  the 
surface  S  =  F{U)  at  each  U .  Let  /J  be  the  components  of  F  and  consider  the  action  of  the  Hessian 
V2F(a,  b )  on  the  vectors  a,  b  €  IRN.  V2F(a,  b )  is  vector-valued  with  components  a  •  V2/-'  b.  Since 
V2/J  is  symmetric,  it  follows  V2F’(a,6)  =  V2F(b, a).  For  U  fixed,  t  €  IR  and  a,  b  €  IR^  equation 
(3.6)  implies 

(3.34)  lk(U  +  ta)  •  VF{U  +  ta)b=  X k(U  +  ta ) lk(U  +  ta)-b. 
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Differentiating  (3.34)  with  respect  to  t  and  setting  t  =  0  in  the  resulting  equation,  we  deduce  the 
identity 

(3.35)  lk  ■  V2F(a,  6)  =  (VAfc  •  a)  (lk  ■  b )  -  (V/*  a)  ■  (VF  -  Xk  I)b , 

which,  in  turn,  yields  the  well-known  identities 

h  •  V2F(rm,  rn)  =  (VAjt  •  rm)  (/*  •  rn)  +  (A*  -  A n)(V/fc  rm  •  rn) 

(3-36)  =  f  (Ajt  -  An)(V/jt  rm  •  rn)  A:  #  n 

\  (VA k  •  rm)  (/fc  •  rfc)  k  =  n 

The  coefficients  f}k,mn  are  related  to  the  second  derivatives  lk  •  V2F(rm,rn)  whenever  k  ^  m  or 
k  ^  n.  There  is  also  the  formula 

(3.37)  (Afc-An)(V4rfc.rn)  =  (VAfe.r„)(/fc.rife), 

The  coefficient  fik,kk  —  V/fc  rk  •  rk  does  not  appear  in  the  above  relations.  To  explain  that 
consider  the  effect  of  renormalizing  the  eigenvectors  on  the  coefficients  f3k,mn  and  especially  to 
0k,kk.  Let  {ffc}  and  {4}  be  a  given  set  af  right  and  left  eigenvectors  and  set  rk  =  rkfk ,  lk  =  sklk 
where  Tk  =  Tk(U)  and  sk  =  sk(U)  renormalizing  factors  with  Tk  >  0,  sk  >  0.  A  simple  computation 
shows  V/fc  =  lk  ®  Vsfc  -(-  skVlk  and  thus 

0k, mn  —  =  Tm^n  [(f*m  *  ^ ^fc)  (4  ’  f"n)  *  f’n] 

(3-38) 

=  Tmrn  [(fm  •  Vsfc)  (4  •  f„)  +  ,mn] 

If  A;  ^  n  the  renormalization  has  no  effect  on  the  sign  of  0k>mn.  If  k  =  n  though  the  renormalization 
of  the  left  eigenvectors  has  an  effect  on  0k>mk  and  can  make  it  to  be  zero. 

In  particular,  if  we  restrict  in  the  neighborhood  of  some  state  U,  we  can  choose  a  renormal¬ 
ization  so  that  the  resulting  eigenvectors  satisfy  simultaneously 

(3.39)  lk‘rk  =  0,  V/jferjt  .ffc  =  0. 

To  this  end,  choose  first  sk  so  that 

(3.40)  (ffc  •  Vsjt)  (4  •  fk)  +  skVlk  fk  •  ft  =  0 . 

(3.40)  is  a  hyperbolic  equation  for  sk.  If  we  give  data  for  sk  in  a  hypersurface  S  transversal  to 
the  vector  field  rk,  the  Cauchy  problem  for  (3.40)  has  locally  a  unique  solution.  If  the  data  are 
positive  then  sk  >  0.  Following  that  rk  is  chosen  so  that  rksklk  •  rk  =  1.  The  resulting  {rk},  {lk} 
have  the  desired  properties. 
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4.  Properties  of  the  functions  (fk  -  Wave  interaction  estimates 
Let  C°(— 00,00)  stand  for  the  space  of  the  continuous,  bounded  functions.  Consider  the  set 


(4.1)  n  =  {vec°(- 00,00):  sup |v(0  —  U-\  <  p} 

£€IR 

and  suppose  that  p  is  sufficiently  small  so  that  the  wave  speeds  A k(V)  are  bounded  and  totally 
separated  for  V  €  ft  : 

(At)  Ajt_  <  A*(V(0)  <  Afc+ 


Ai-  <  Ax (y(0)  <  Ai+  <  A2_  <  A2(F(0)  <  A2+  <  ... 

( A<i ) 

<  A(at_i)_  <  Ajv-i(nO)  ^  A(at_i)+  <  A <  Ajv(V(0)  -  Atv+  • 

Consider  the  linearized  equation 

(4.2)  £¥>'fc  +  (e-  Ajfc(F(O))^  =  0. 


The  fundamental  solution  of  (4.2)  may  be  written  in  the  form 


(4.3) 

<Pk  = 

where 

(4-4) 

II  II 

<0 

e~isk  _  1  -if(  s-\k(V(s)) 

fZ*-iet(c)dC  ~^e 


ds 


I! 


d(  ■ 


Recall  that  gk  has  the  form  of  a  potential- well  function  ( c.f .  Figure  1)  and  that  pk  is  selected 
as  a  point  where  gk  achieves  its  global  minimum.  As  a  result  pk  satisfies  A k-  <  pk  <  A^+, 
\k{V{pk))  -  Pk  and 


(4-5)  9k(O>9k(Pk)  =  0,  £  €  IR. 

The  intent  of  this  section  is  to  establish  various  estimates  on  the  functions  <fk  and  integrals  involving 
them,  that  are  needed  in  the  forthcoming  constructions. 

We  begin  with  a  careful  analysis  of  the  behavior  of  (fk  in  the  limit  £  — *  0.  Given  a  positive 
function  h(e),  we  use  the  customary  notation  f(e)  =  0(h(e))  as  e  — >  0  to  mean  there  are  constants 
£0  sufficiently  small  and  C  such  that  |/(e)|  <  C  h(e)  for  0  <  e  <  £o- 
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Lemma  4.1.  Suppose  the  wave  speed  A*(V)  satisfies  (Ai). 


(i)  7/  4  =  A/t+  -  A*,. 

>  0,  then  as  e  — ►  0  : 

(4.6) 

(4.7) 

o  <  <pk(0  <  0(1)  y  ,  forfelR, 

and 

Vfc(0  <  0(1)  ^  ,  for  £  <  A* _  , 

(4.8) 

£ 

¥>*(0  <  0(1)  y  e-*«"A*+>2  ,  for  f  >  Afc+  . 

(ii)  7/  =  Afc+  -  Afc_ 

.  =  0,  then 

(4.9) 

he  =  V2tt£  ,  <Pk(0=  /K - e  2cU  A*_)  • 

V27T6 

Proof.  Assume  first  that  <4  >  0.  Performing  the  change  of  variable  (  =  pk  +  \/^  V  in  the  integral 

(4.4),  we  obtain 

roo  roo 

Ike  =  /  e" * Sfc(c)  dC  =  \/£  / 

(4.10) 

7  —  00  7  —  00 

roo  _ij'k+**,_Xi(v(.))d 

=  y/e  1  e  Jp*  dt]. 

J — oo 

Using  again  the  change  of  variable  s  =  pk  +  y/er  and  (4.5),  we  have 
1  1 

-  9k(pk  +  y/er])=  ~  *  -  (5))  ds 

£  £  Jpk 

(4'n)  =  Jo  {r  "  [Xk{v{pk  +  ^r))  -  A*(v(/>*))]  }  dT 

>  0  ,  for  Tj  £  IR . 

We  remark  that  for  r;  >  0 

(4-12)  £  T  -  -L  [Afc(V(pfc  +  y/i  r))  -  \k(V(pk))  ]  dr  <  £  +  -L  (Afc+  -  Afc_)  7 7 

while  for  rj  <  0 

(413)  J*  T~7!  ^k(y^k  +  v£r))  -  A»(V(«))]  dr  <  £  -  -L(A1+  -  Xk-)ij. 
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Therefore  (4.10),  (4.11),  (4.12)  and  (4.13)  provide  the  estimate 

he  -  Vs  J  exp  |  -  J  t  -  —jz  [A/t(F (pk  +  \fi r))  -  A k(V(Pk))  ]  dr j  dr, 
+  y/e  J  exp  |  -  J  t  -  [Afc(V(pfc  +  V^T))  ~  Afc(V(pfc))]  dr|  dr, 

(4'14)  >  yfe f  dV+y/Ie%  I™  dr, 

7-oo  Jo 

The  asymptotic  behavior  of  the  last  integrals  can  be  evaluated  by  using  the  limits 


r°°  —  £-  is¬ 
le  *  dc 

lim  — - 5 —  =  hm 


(4.15) 


—e 


*2 

2 


=  1 


1  -  £- 
~e  2 

X 


x* 


e  2 


(4.16) 


Jx  e-*Td(  e~4 

Urn  — - =  lim  - - p - -r  =  1 

*->-<*>  _Ie-V  x-f-oo  4-  e“V 

and  yields  for  small  e 
(4.17) 


.0(1)  dk 

Next,  observe  that  for  £  >  Xk+  >  Pk 

K 


0(1)  dk 


(4.18) 


while  for  £  <  \k-  <  pk 


9k(0=  \  s  -  Xk(V(s))ds  + gk(\k+) 
hk+ 

>  f  (s-Aw)*=i«-At+)2, 
J  Al+  1 


(4.19) 


9k(0  =  t  s-\k(V(s))ds  +  gk(\k-) 

Jx  k- 


>  -  J  "  (s-Xk-)ds=  i(e-  Afc_)2. 


Therefore,  (4.4)  and  (4.5)  imply 

r*t- 


(4.20) 


/Ak_  f  oo 

e-±(C-\k_f  d£  +  dk+  I  e-^(C-A*+)  d( 

•OO  J^k  + 

/oo 

e_i’7  dr,  =  dk  +  V2tTe 

-OO 
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which  together  with  (4.17)  complete  the  proof  of  (4.6). 

Estimates  (4.7)  and  (4.8)  follow  from 

w,  ,  ... 

(4.21)  MO  =  —r - <  0(1) -J  e- , 

*ke  £ 

a  consequence  of  (4.3)  and  (4.6),  in  conjunction  with  (4.5),  (4.18)  and  (4.19).  Finally,  if  dk  =  0  then 
Afc(F)  remains  constant,  say  \k-i  and  (4.9)  follows  from  (4.3)  and  (4.4)  via  a  direct  calculation.  | 

Remark.  To  expand  on  the  implications  of  the  lemma,  suppose  we  are  given  a  family  of  functions 
{Ue}e> o  C  ft  and  that  for  each  Uc  we  define  the  corresponding  <pke  solution  of  (3.24).  Then  (4.8) 
implies  that  <pke  — 1 ►  0  as  e  — ►  0  uniformly  on  any  interval  of  the  form  (-00,0*]  U  [5*,  00)  with 
a.k  <  Afc_  <  A k+  <  bk.  The  family  {<Pkc}c>o  is  uniformly  bounded  in  L 1  and  thus  there  exists  a 
subsequence  <pkCn  with  £„  -*•  0  and  a  finite  Borel  measure  <f>k  with  supp  <f>k  c  [A|<_,  Ak+]  such  that 
<Pkc„  — 1 ■  4>k  weak-*  in  measures.  For  the  single  conservation  law  or  the  equations  of  isothermal 
elasticity,  objects  similar  to  <f>k  yield  the  same  structure  for  the  Riemann  problem  solution  as  that 
obtained  by  the  Liu  shock  admissibility  criterion.  ( cf  Tzavaras  [TZ2]). 

Our  next  task  is  to  study  certain  integrals  involving  ipm  and  <pn  that  calculate  the  effect  of 
interactions  between  elementary  waves.  It  is  convenient  to  introduce  the  notation: 
dk  =  length  of  the  interval  [Afc_,Ajt+] 

Cfc  =  middle  point  of  the  interval  [A/-_ ,  Ajt+] 

d(€,Xk)  =  distance  between  the  point  £  and  the  interval  [Afc_,Ajt+] 

Dmn  =  d( Am,An)  =  distance  between  the  intervals  [Am_,Am+]  and  [An_,An+] 

Because  of  (/1 2)  Dmn  >  0;  also,  we  may  assume  without  loss  of  generality  that  dk  >  0  by  replacing 
(4.9)  with  the  weaker  estimates  (4.6  -  4.8).  Lemma  4.1  indicates  that  <pk  has  the  form  shown 
in  Figure  2.  The  behavior  of  <fk  is  uncontrolled  in  the  interval  [Afc_,Ajt+],  where  the  wave  speed 
A*(F)  takes  values,  but  its  amplitude  is  of  at  most  order  0(4).  For  f  £  [Afc_,  Afc+],  <pk  decays  like 
0(j  e~^d^,xO  ).  it  is  expedient  to  fix  points  a*,  bk,  k  =  1,  ... ,  N,  such  that 
<  Aj_  <  Ai+  <  b\  <  a2  <  A2-  <  A2+  <  62  <  ••• 

(4.22) 

<  tiN-i  <  A(^r_!)_  <  A t tv — 1 )+  ^  bpr- 1  <  ® n  <  A n-  <  Ayv+  <  hjv  , 
and  introduce  the  notation 

(4.23)  sk(e)  =  max  <pk(£) ,  ak  =  ^  minfla*  -  Afc_|2,  \bk  -  Xk+\2}  . 

Then  (4.8)  implies 

(4-24)  MO  <  *k(e)  <dkO(-ee~^),  [ak,bk] , 
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1  **  k 

and  Sk(s)  serves  as  a  global  bound  outside  the  main  support  of  the  wave.  The  function  =  -  e  *  , 
describing  the  decay  rate  as  $  — »  0,  has  the  behavior:  hk  increases  from  0  to  its  maximum  value, 
achieved  at  £  =  a*  and  equal  to  1  / e  a*,  and  then  decreases  down  to  0  as  e  — ►  oo. 


FIGURE  2. 


Lemma  4.2.  Suppose  the  wave  speeds  Afc(V')  satisfy  (.A2)  and 


(.4.3)  (1  +  \/3)  ( dm  +  dk)  <  d(Afc,  Am)  —  Dkm 

for  U  6  fi.  Then  there  exist  constants  a  km  >  0  depending  on  dk,  dm,  Dkm  but  independent  of  e, 


V  such  that 


(4.25)  |e_*“  /  e‘5*«Vm(C)<*C|  <{,?*”*  £^m+0^e  ^  *’ 

1  JCk  1  \\t-Ck\<Pk,  m  =  k. 

Proof.  When  m  =  k,  (4.25)  follows  from  a  direct  calculation.  So  suppose  that  m  ^  k.  Using  the 

notation  A k  —  Afc(V(f))  and  (4.3)  we  obtain  the  chain  of  identities 


’is*  e*i 

J  Ci- 


>Vm(Od(  =  e  tJn 


I  s- A*  ds  ft  }  Jc  s-\k 
«  JfiL  I  0* 


f  .  pf  ~~  \  f  a  Am  ds 

*  1  P  .9 — At  ds  e 


(4.26) 


j  *  ^>m _  ft  *~Xk  ds  ft  s~Xmds , 

I  me  Jci 


J  Ct 


In  view  of  (A2 ),  we  have 


(4.27) 


/  e‘f  Am_A*dsdC  <|±/  f  fc| 

J cjt  I  */  cjt  u(Am,A*) 

>  I .  'J1- 1 ‘ 

Dmk  I  Jck  J  £ 


■Dmfc  |  7ct  ' 

.  £  A  ,  -i/£  Am-A*dA 

sI^l1+e  *  j 


Am  -  Afcds 
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Combining  (4.26)  with  (4.27)  and  using  (4.3),  we  arrive  at  the  estimate 


The  goal  is  to  show  that  under  (A3)  the  term  in  parentheses  decays  as  e 
observe  that 


0.  To  this  end, 


(4.29) 


s-  Xmds<  --D 


2 

km 


(4.30) 

It  suffices  to  show 

(4.31) 


r 

J  Pk 


Am  —  Afc  ds  dk  ( Dkm.  "I"  dk  +  dm) 


— ojkm  :=  —-D\m  +  (dk  +  dm )  Dkm  +  (dk  +  dm )2  <  0 


Since  the  roots  of  the  quadratic  —  jx2  +  x  +  1  are  1  ±  \/3,  hypothesis  (A3)  implies  the  inequality 

(4.31)  and  thus  there  exists  a  positive  constant  atkm  such  that 

-  j  fPk  s-\m  /c*  -A*  ds 

t>  1  JPm  1  JPk 


(4.32) 


<  0(e  *  ) . 


The  proof  of  the  lemma  follows  from  (4.28),  (4.32)  and  (4.6).  | 

Our  next  objective  is  to  use  the  facts  that  each  <pk  is  essentially  supported  on  the  interval 
[A*,-,  Afc+]  and  that  such  intervals  are  distinct  in  order  to  estimate  the  integrals 

(4.33)  Fk,mn( 0  =  /*  «‘*(CW0  MO d(. 

Jck 

We  begin  with: 

Lemma  4.3.  Suppose  that  Ajt(Vr),  k  =  1,  ... ,  N ,  satisfy  (A2)  and  (A3).  Then 
(i)  for  m  =  1,  ... ,  N, 

(4.34)  |  e~‘9k  f  el9k  <pm(pkdC\<<Pk, 

Jck 
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(ii)  for  m,  n  =  1,  ... ,  N,  with  m  ^  n,  m  k  and  n  ^  k, 


(4.35) 


\e  ‘ 9k  f  e^3k  <pm  (pnd(  |  <  <pn  +  ~  rT ~~  <fr 

J  CJ. 


Dkn 


Dkr, 


r  /  \  s\/  ***"  v  dkdn  t  \  r\t  —  *k,m.  \  dkdm  ■< 

+  [sm(e)0(e  «  )— —  +  sn(e)0(e  «  )-~ — j  <fk 

L'kn  '■ 


Dkm 


Proof.  First  we  show  (i).  Since 


(4.36) 


1  ft  1  e-*3*  ft 

Fk,mk  =  e  ‘9k  /  e‘9k(pm— — d(  =  <pk  <pmd<;, 

J Cjt  Cjb 


it  follows  |Ffc,mfc|  <  and  (4.34)  is  proved.  Observe  next  that  because  of  (A^) 


(4.37) 


<Prn  <Pn  <  sm{s)  <Pn  +  sn(e)  ipm ,  form/n,  (el. 


Using  (4.25)  with  m  ±  k  and  n  ^  k,  we  obtain 

r€ 


(4.38) 


\Fk,mn\  <  «m(e)  |e  ‘at/  e<3*  <fn  d(  \  +  s„(f)  |e  ‘3*  f  e «  3fc  V7m  I 

Jcjt  •'c* 


< 


^m(^) 


(e^n  +  0(e  ~J^L)dkdn(pk)  +  (ey>m  +  0(e  _^IL)  dkdm  <pk) , 


which  in  turn  yields  (4.35).  | 

It  remains  to  estimate  the  integrals  Fk<mm  with  m  /  fc,  that  account  for  the  effect  of  self¬ 
interactions.  Using  (4.3),  we  write  Fk<rnm  in  the  form 

r€ 


Fk,  r 


—  e  t3k  f  e\gk  ^,2 

Jck 


<Pmd( 


(4.39) 


fC  *-\kds  ft  j  f(  s-\kds  e  *  m 

z=e  ,J*k  /  ^ - d( 

Jck 


'Ck 


1 2 


e-?  JL  '-X”  da  ft  J  f<  ,_Xk  d3  _2  f<  a.Xti 

~  p  J  e  e 

■Lme  J  Ck 

2  /"£  ~e  ft  s~^km  d.3  .  . 

=  <Pm  e  cJ(  d( , 

«/  Ci 


da 


dC 


where  we  have  set 


(4.40)  A*m(U)  =  2  Xm(U)  -  A *(17)  =  Am(U)  +  (Am(ff)  -  Afc({/)) . 

Note  that  the  ordering  goes  A*(U)  <  Am(U)  <  Akm(U)  when  k  <  m  and  A km(U)  <  A m(U)  < 
Xk(U)  when  k  >  m.  In  order  to  estimate  FkiTnm,  it  is  necessary  to  study  the  ranges  of  the  wave 
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speeds  Afc(F)  and  Am(V)  relative  to  the  range  of  the  composite  speed  Afcm(V),  for  V  €  ft,  and 
to  impose  conditions  that  guarantee  non-resonance  between  the  wave  speeds  and  the  composite 
speed.  Note  that  Afcm(Vr)  is  bounded  by 

(4.41)  A*m_  <  Afcm(V(f))  <  Afcm+ 

where  the  constants  A*m_,  A *m+  and  dfcm,  the  length  of  the  range  of  Afcm(F),  depend  only  on  /i. 
We  introduce  the  notation 

<*((,  A*km)  —  distance  between  the  point  f  and  the  interval  [Ajtm_,  Afcm+], 
d(Xm,Akm)  =  distance  between  the  intervals  [Am_,Am+],  [Ajtm_,Afc  m+]j 
and  impose  the  strengthened  version  of  Hypothesis  (A3) 

(A4)  7  (dTO  -J-  dk')  —  7  [(Am4-  —  Am_ )  -4-  {Xk+  ~  A^_)]  <  d(A^,  Am)  =  Dkm  • 

It  is  easy  to  calculate  dkm  =  Akm+  -  Akm-  =  2 dm  +  dk ,  d(Am,Afcm)  =  Dkm  ~  dm,  and  to  note 
that  the  ranges  of  A^(F),  Am(F)  and  A*:m(Vr)  are  separated  for  V  €  D  (see  Figure  3). 

A k(V)  Am(V)  Akm(V)  =  Am(F)  +  (Am(V)  -  Xk(V)) 

— ( - ) - 4 - ) - ( - )- 

Q>kXk—  X  k+bk  dm  Am_  Xm+bm  dkm  A  km—  Akm+  ^fem 

FIGURE  3. 

The  ranges  of  Ajt(F),  Am(V)  and  Ajtm(F)  for  m  >  k. 

Since  the  lengths  dk  are  of  order  O(fi)  while  the  distances  Dkm  are  of  order  0(1)  as  /z  — *  0, 
hypotheses  (A3)  and  (A4)  are  not  particularly  restrictive  for  solutions  of  small  oscillation.  (A3— A4) 
are  imposed  for  all  k,  m  =  1,  ... ,  N,  and  points  a*m,  are  selected  (near  the  support  of  Afcm(U)) 
so  that,  upon  rearranging  am,  bm  if  necessary, 

(4.42a)  Ok  ^  Xk-  S  Xk+  K.  bk  dm  Xm—  —  Am.f.  <C  bm  dkm  ^  A  km—  —  A  km+  ^  bkm  > 
when  k  <  m,  and 

(4.426)  dkm  A  km—  ^  A  km+  ^  bkm  ^  dm  <  Xm—  ^  Am_|_  <  6m  <  Ok  <  A  k—  ^  A/.4.  <  bk  , 

when  k  >  m.  Such  choices  are  clearly  possible.  The  points  a*,  bk  are  now  fixed,  while  the  points 
akmi  bkm  will  be  selected  subject  to  (4.41)  in  the  course  of  proving: 
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Lemma  4.4.  Suppose  that  Afc(F),  Am(F)  satisfy  (A2  —  A 4).  There  exist  choices  of  akm,  bkm  and 
constants  a*m,  (3km  >  0,  depending  on  dk,  dm,  Dkm  but  not  on  e,  such  that 

(a)  if  k  <  to,  then 


(4.43a)  |e  >Sk  f  e^3k  d£|  < 
(b)  if  k  >  m,  then 

(4.436)  f*  e±9k  cp2md<;\  < 

J  Ck 


d(akm\x^)  £  Tm  +  0(ee  ‘  )  W  >  ^  °*”1’ 

dkmdm0{\e  ^a')y>m  £  >  akmi 


,x  _2«km 


dkm  dmO(\  e  ^°')y>rn  Z<bkm, 

dib^M-)  £  Vi  +  °(*e  t>hm, 


Proof.  Let  k  <  m  and  proceed  to  prove  (a).  The  ranges  of  Afc(F),  Am(F)  and  A fcm(F)  for 
V  €  D  are  as  in  Figure  3,  and  akm  is  any  point  compatible  with  (4.42).  Let  pkm  be  a  point 
where  the  function  J^s  —  Akm{V(s))ds  achieves  its  global  minimum.  Then  A km(V(pkm))  =  Pkm , 
Akm—  —  Pkm  —  A  km+  and 

(4.44)  Gkm(0  =  f  S  -  A  km{V{s))  ds>  0 ,  forwent. 

J  Pkm 

Consider  first  the  region  f  <  akm  <  A  km--  In  this  region  Fk,mm  in  (4.39)  is  decomposed  into 
the  integrals 


„  2  c  f(  3~Akm  ds  ( 

(4.45)  Fh,mm  =  <Pme  J '**  (  J 


-i/(  3-Akmds  [Ck  -j  A*„ 

‘  J»km  dC  —  /  e  Jpkm 


■  /Ct  -i  /'  *-A kmds  \ 

dC  ~  e  Jpkm  dC,  J 


The  first  integral  is  dominant  when  £  >  ck  and  the  second  is  dominant  when  £  <  ck.  Since 
c<£<  akm  <  Afcm_,  the  first  integral  is  estimated  by 


L 


^  -j  f<  »-Atm  ds 

e  Jpkm 


r€ 


it.  a-A*m  rfa  ( Afcm  C 


|>  ^  /  —7  /  O®  / 

d(<  e  J'*»*  ( 


(4.46) 


00 

€ 


Akm  £ 


)d( 


ds 


_£  ^  *~A*1 

€  ;rA  -T - dl--/  s  —  Akm  ds 

d(£,Akm)  V  £Jpkr 


s~  Afcm rfs) 


-7  /* 

3  * 


ds 


Afcm) 

In  a  similar  fashion  the  second  integral  is  estimated  by 


(4.47) 


/_ 


C*  _1 


if*  -A*. 


ds 


< 


A/*m) 


j  fCk  S-A*m  <*S 


£  -7  fCk  s-Akm  ds 

<  - - - - -  e  JPJfcm 
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Let  D\  =  d(afcm,  Afcm),  D2  =  d(Xk,  Afcm)  and  combine  (4.45),  (4.46),  (4.47),  and  (4.3)  to  obtain 


lr>  -  i  (  £  £  if*  *-A*m  da\ 

\Fk,mm\<<Pm  \j^  +  “  J 


(4.48) 


-1  r(  5—  Am  ds 

e  2  ee  A 

=  &f<Pm+D~2  /L- 


}  f *k  a—Akm  ds  A  f4  s-Atm  ds 

e  Jc*  e  Jp't 


d1  +  d2  n 


£  1  /  -f  fP*  «-Am  ds+}  ffik  *-Akm  ds 

_  _  *  Jpm  C  Jet 


'  Am 


)  he  <Pk  ■ 


me 


It  suffices  to  show  that  the  term  in  parentheses  decays  as  e  — *•  0.  Using  the  estimations 

rPk 
'Pm 


Am  ds  <  —D\ 


(4.49) 


fPk 

-2  S- 
Jpm 
fPk 

I  $  ~  A- km  ds  <  (A km+  A/;_)4  —  2 [D^m  4*  (dm  +  4)](4i  +  4)  , 

•'Ct 


together  with  the  fact  that  (A3)  implies  that  (4.31)  is  satisfied,  we  conclude  that 

_2  rfik  s-\m  ds+  \  ffik  s-A km  ds 

0  J  Pm.  C  Jck 


(4.50) 


<0(e"=fflt). 


In  conjunction  with  (4.48)  and  (4.6),  (4.50)  shows  (4.43)  for  f  <  dkm,  k  <  m. 

Consider  now  the  region  £  >  a^m.  An  argument  similar  to  the  one  leading  to  (4.20)  shows 

that 

(4.51)  |  f(e'ih-A‘-d‘d(\  <  f 

J  Cjfe  ./  —  c 

Therefore,  (4.39)  and  (4.40)  give 


— -J  f(  3-Alcr, 

o  c  Jp 


ds 


<*C  ~  dfern  +  V 2  7T  £  . 


I^.mrnl  <  0(1)  4m  ^  C  ‘  *  A*~  * 

/4  (-o'!  /0/i\j  e  *  “  m  t  J*  s-Xmds  -jfe  \m-Akds 

(4.52)  =  0(1)  4m  - ~T - etJ'km  e  Jnm 


<  0(1)  4m  fm  ^ 


dm  —1  f  km  s— Am  ds  J  J  Am—A*  ds 
-  e  C  JPm  e  Jfikm 


The  goal  is  to  choose  a*™,  so  that  the  term  in  parentheses  decays  as  £  — ►  0,  for  any  £  >  a^m.  Since 
^  [Am^,Am+],  the  first  term  decays  as  e  — ►  0  and  its  decay  rate  can  be  estimated  by  noting 

that 

r  Ai-m_ 


(4.53) 


rpkm  [Akm- 

—  I  s  -  \mds  <  -  I  s-  Am+  ds 

J  Pm  j \m  + 


-^(Afcm-  -  Am+)2  =  ~\{Dkm  ~  dmf  . 
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Since  Xm(U)  >  A^.( U),  the  second  term  decays  for  £  >  pkm  but  grows  for  £  <  Pkm ■  The  fastest 
growth  occurs  for  £  =  Ofcm  and  the  growth  rate  is  estimated  by 


(4.54) 


Am  —  A kds  <  (Am+  —  A k—){Pkm  akm) 


<  ( Dkm  +  dm.  +  dk)  (d(a,kmi  A-km)  +  2  dm  +  dfc)  . 

It  suffices  to  give  conditions  on  dk,  dm,  Dkm  and  to  choose  akm  so  that 


(4.55)  ~Pkm  '■=  — “ (Dkm  —  dm )2  +  ( Dkm  +  dm  +  dk )  (d(afcm,  Afcm)  +  2  dm  +  dk)  <  0  . 


For  example,  if  we  choose  =  Akm-  —  dk  and  we  require  that 


(4.56)  4  ( Dkm  +  dm  +  dk)  ( dm  +  dk)  <  [Dkm  —  {dm  +  ^fc)]  > 


then  (4.55)  is  satisfied.  By  solving  the  inequality  y 2  -  6 xy  —  3x2  >  0,  we  see  that  (At)  implies 

(4.56) .  Therefore  (4.52)  yields  the  estimate 

X  ft 

(4.57)  |^fc,mm|  <  dmO{—  C  c  )  tym 

for  £  >  ajtmj  k  <  m,  and  completes  the  proof  of  part  (a).  The  proof  of  part  (b)  is  similar.  | 

Lemmas  4.3  and  4.4  provide  estimates  on  the  integrals  Fk}mn >  which  calculate  the  effect  of 
diffusion  induced  wave  interactions.  The  estimates  are  consequences  of  the  separation  hypotheses 
( A2  —  At)  on  the  wave  speeds.  Obviously  (A*)  is  the  strongest  hypothesis  and  implies  the  rest.  In 
the  sequel  we  make  use  of  the  following  implication  of  (4.34),  (4.35),  (4.43)  and  (4.7). 

Corollary  4.5.  Suppose  that  Ajt(F)  satisfy  (At),  for  ky  m,  n  =  1,  ... ,  N.  There  is  €q  >  0  and  a 
constant  C,  depending  on  <4,  Dkm,  D^n  but  not  on  £,  such  that 

ft  N 

(4.58)  \Fk,mn\  =  \e~‘9k  /  e**  <<?$>; 

Jpk  j= 1 

for  k,  m,  n  —  1,  ... ,  N  and  0  <  e  <  £q. 


Remark.  It  is  instructive  to  identify  which  of  the  integrals  Ffc,mn  have  nonzero  contributions  in 
the  limit  e  — ►  0.  In  view  of  (4.36)  and  (4.3),  the  terms  Fk,mk  and  Fk,km  Fk,kk  have  nonzero  limiting 
contributions  supported  on  the  fc-th  wave  speed.  On  the  other  hand,  (4.35)  and  (4.7)  imply  that 
Fk,mn  -  0  as  £  — ►  0  when  m  ±  n,  m  ±  k  and  n  ±  k,  what  suggests  that  diffusion  induced 
interactions  of  two  distinct  families  have  no  contribution  as  £  —*■  0  on  a  third  family.  (Recall  that 
we  are  dealing  with  Riemann  data  solutions).  By  contrast,  (4.43)  suggests  that  the  terms  Fk,mm , 
m  ^  k,  accounting  for  the  effect  of  self-interactions  on  another  family,  have  a  nonzero  contribution 
in  the  £  — *  0  limit  which  is  supported  on  the  m-th  wave  speed. 
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5.  Validation  of  the  asymptotic  expansion 
The  objective  of  this  section  is  to  solve  the  problem 

«*»+ [e  -  A*(v(«)} 

N  N 

(5.1)  =  £  E  E  [VAfc(F(f))  rm(V( 0)  •  rn(V(0))  (rm  <fim  +  0m)  (r„  +  6n) 

m=l  n=l 

^fc(cfc)  =  0 

where  V  e  O,  defined  in  (4.1),  and  r  =  (ri,  ... ,  tn)  is  a  vector-parameter  in  IR£  The  aim  is  to 
construct  solutions  9k( •  ;r)  that  are  of  order  0(|r|2)  in  the  wave  strength  |r|  =  1^1  +  ...  -f  |r/yj  as 
|r|  -+■  0.  This  would  validate  the  asymptotic  expansion  (3.29). 

Throughout  the  section  we  use  the  notation 

A*  =  Afc(F(0), 

(5.2) 

Pk,mn  =  Pk,mn(V(0)  =  VA k(V(0)  rm(V(0)  ■  rn(V(0) 

and  assume  that  //  is  small  so  that  the  hypotheses  (A:  -  A4)  on  the  wave  speeds  are  fulfilled  for 
V  €  ft.  Moreover, 


(5.3) 


\Pk,mn\  <  B 


with  B  depending  only  on  fi.  Recall  that  gk  is  defined  in  (4.4)  and  that  ck  is  the  middle  point  of 
the  interval  [Afc_,Afc+].  Using  the  variation  of  parameters  formula,  (5.1)  is  expressed  as  a  system 
of  integral  equations 

ti  N 

(5.4)  9k(0  =  e~‘9k  /  £  Pk,mn(V(Q)  (rm  <pm( 0  +  MO)  fc  MO  +  *n(0)  d(  ■ 

^ Cfc  m,n=l 

Our  strategy  is  to  formulate  (5.4)  as  a  fixed  point  problem,  and  to  use  the  uniform  contraction 
principle  in  order  to  construct  solutions  9k{- ;  r),  k  =  1,  ... ,  N. 

Let  Co(IR)  stand  for  the  continuous  functions  that  decay  to  zero  as  |£|  — ►  oo  and  define 

(5.5)  E={x  =  (Xu-,XN)e[C0m]N  =  SUP  <  00 ’  i  =  h-,N 

1  E;=1  Viifl 

E  with  the  weighted  sup-norm 


(5.6) 
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with  weight  <Pi  >  0  is  a  Banach  space.  Let  Bg  =  {r  6  Dlw  :  |r|  <  S  }  and  set 
'  N 

(5.7)  F={XeE:  l*i(f)l  <  A  |r|s  ?€lR.j  =  l . N  }  , 

«=1 

where  r  €  B$  and  A  is  a  constant  to  be  determined  later.  F  is  a  closed  bounded  subset  of  E  in 
the  weighted  norm  ||  •  ||.  Define  the  map  T  that  takes  V  €  O,  t  €  Bg,  X  €  F  to  the  vector-valued 
function  T(x)  with  components 

if 

(5*8)  Tk(x)  =  €  %‘9k  I  €~3k  ^  ^  fik,mn  (^m  V’m  4"  Xm)  (Tn  ‘Pn  4"  Xn)  ) 

•^c*  m, n=l 

fe  =  1,  ... ,  N.  The  map  T  has  the  properties  : 

Proposition  5.1.  There  exist  positive  constants  A  and  So  such  that  for  6  <  So  : 

(i )  T  :  Q  X  Bg  x  F  —*  F  is  well  defined. 

(ii)  There  exists  a,  0  <  a  <  1,  such  that 

(5.9)  ||  T(V,  r,  X)  -  T(V,  r,  x)  ||  <  <*  ||X  -  *11  ,  for  X,  X  €  F , 

and  for  any  V  €  D,  r  €  Bg.  Therefore  T(V,t,  •) :  F  — >  F  is  a  uniform  contraction. 

(iii)  There  exists  a  positive  constant  C,  depending  on  p  but  independent  of  6,  such  that 


(5.10) 


||  T(V,  r,  X)  -  T(V,  s,  x)  ||  <  C  6  \t  -  s\  ,  for  t,s€Bs, 


and  for  any  V  £  ft,  x  £  F- 

Proof.  In  the  forthcoming  estimates  C,  C'  and  C"  stand  for  generic  constants  that  can  be 
estimated  in  terms  of  B,  the  dimension  of  the  system  N,  and  the  constant  in  the  estimate  (4.58). 
As  a  result  such  constants  ultimately  depend  on  p  in  (4.1),  but  are  independent  of  6.  We  proceed 
to  establish  (i).  Let  V  €  ft,  r  €  Bg  and  x  €  F  be  fixed.  Then  (5.8),  (5.7)  and  (4.58)  imply 

,  I  /■€  N 

|T*(x)|  <  e“«'* 


<  5e  c9k 


(5.11) 


/  e^9k  53  l^.mnl  (|rm|^m  +  |Xm|)  (|T-n|^n+  |Xn|)^C 
•'c*  m,  n=l 

/  e‘3t  53  {\Trn\<Pm  + A\T\2YjVi)  {\Tn\  <Pn  +  A\r\2  Y  <Pj)  dC 

c*  m  ,n=  1  »‘=1  j— 1 

N  ri 

<  C  |r|2  (1  +  2AS  +  A?  S'1)  Y  \e~‘9k  ei»<pm<pnd(  | 

m,n=  1  Jck 
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Comparing  the  outcome  with  (5.7)  we  see  that  if 


(5.12) 


C(l  +  A6f  <  A, 


then  T(  V,  r,  x)  €  F  and  (i)  is  established. 

Next,  we  examine  (ii).  Let  V  €  fi,  r  €  B$  be  fixed,  and  consider  x>  X  €  F.  Then 


N 


,  .  Tfc(x)  —  Tfe(x)  —  e  ‘9k  f  ec9k  ^  ]  0k,mn\'rm(Pm{Xn  Xn)  +  Tn  tfin  (Xm  Xm) 

V5'13)  Jck  m,n=l 

+  (XmXn-XmXn)]dC 


Using  (5.6),  (5.7),  (5.11)  and  (4.58)  we  obtain 

r€  N 


(5.14) 


I  n 

|Tfc(x)  —  Tfc(^)|  <  e  I  e,9k  ^  ^  |/?fc,mn|  [l^Vnl  |Xn  —  Xn  |  +  |Ui|  V’n  |Xm  —  Xm| 

1  Jc*  m,n=  1 

+  IXm  |  |Xn  -  Xn  |  +  |Xn|  |Xm  “  Xml]  <*C 
I  fi  N  N 

<Be-i9k\  e^9k  ^kmkmllx-xll  ^<Pi 

I  ^ ck  m,n— 1  «'=! 

+  2  A  |r|2  (^2  ipj)  ||x  —  xll  53v><]  dC 

i= l  i=i 

<  C1  (6  +  A62)  (  /  e^9k(PmlPn  d(  |  )  llx  xll 

'  m,n=X  */c*  / 


AT 


<c'«(i  +  /i«)(I>;)llx 

j=l 

which,  on  account  of  (5.6),  in  turn  implies 


(5.15) 


l|T(X)  -  T(x)||  <C' 6  {\  +  AS)  ||x  -  xll  • 


Therefore  T  will  be  a  uniform  contraction  on  F,  provided  that 


(5.16)  C"tf(l  +  ;M)=:a<  1. 

Note  that  (5.12)  and  (5.16)  can  be  simultaneously  satisfied  for  many  choices  of  A  and  6.  In 
the  sequel,  we  fix  A  =  AC  and  6  <  So  =  min{^,  ^c7).  For  these  choices,  1  +  AS  <  2,  both  (5.12) 
and  (5.16)  are  fulfilled,  and  the  proof  of  (i)  and  (ii)  is  completed. 
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Last,  we  turn  to  (iii).  Let  V  £  ft,  X  €  F  be  fixed  and  consider  t,  s  £  Bs-  Then  (5.8)  yields 
(upon  suppressing  the  x  and  V  dependence) 

fl  N 

Tjfc(r)  —  Jfc(s)  =  e  *  3k  I  ec3k  ^  '  A,mn  [(rm7,n  "  ^m^n) 

(5-17)  Jck  mn=1 

+  (Tto  —  sm )  V’m  Xn  T  (Tn  ~  sn)  Pn  Xm  ]  • 

Using  (5.6),  (5.7),  (4.58)  and  (5.16)  we  deduce 

N 


\Tk(r)-Tk(s)\<e-i» 


ft  ( 

/  e‘3k  \Pk,mn\  (  [Irm  ~  •sm|  |rn|  +  \Tn  ~  5n|  |-Sm|]  Pm  Pn 

Jc*  m,n= 1  ' 

T  l^m  &m |  Pm  |Xn|  T  I^Vi  Pn  |Xm|) 


(5.18) 


l  ri  N 

/ 

/  E 

U[|rm-sm|  +  | 

1  Jck  m,n= 1 

\ 

N 

4~  A6^  (^  ^  Pj)  [l^m  ~  sm  1  Pm  ~H  [ Tn  ^nlV’n]) 
3= 1  7 

N  fi 

<  C" S  (1  +  AS)  |r  -  s|  ^2  |e~'5*  /  ^l9kPmPn  dC| 

_ 1  Jck 


dC 


m,n= 1 


N 


<C"6\r-s\ 

3= i 


and,  by  virtue  of  (5.6), 
(5.19) 


||^(r )  _  T(s)\\  <  C"  ^  |t  -  s| 


which  completes  the  proof  of  (iii).  | 

The  properties  of  the  map  T  are  useful  for  both  solving  problem  (5.1)  and  for  establishing 
properties  of  the  constructed  solution  0  =  (0i,  ... ,  On)- 

Corollary  5.2.  Let  A  and  6  be  as  in  Proposition  5.1.  Given  F  6  fi,  r  G  there  exists  a  unique 
solution  0(-;r)  of  (5.1)  in  the  class  of  functions  satisfying 

N 

(5.20)  |0*(-;r)|  <  A|r|2  |r|<£,  k  =  l,...,N. 

3  =  1 

Moreover,  there  exists  a  constant  C  independent  of  S  such  that  6(-  ;r)  satisfies 

N 

(5.21)  \h(-;T)-ek(--,s)\<C6\T-s\Y,P3,  forr,seBs. 

3  =  1 
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Proof.  For  each  fixed  V  G  ft,  r  G  B&,  the  map  T(V, r,  •)  :  F  —*•  F  is  a  contraction  with  a  uniform 
contraction  constant  a  <  1.  The  first  part  of  the  lemma  is  a  direct  consequence  of  the  contraction 
mapping  theorem. 

The  fixed  point  0  depends  parametrically  in  V  and  r.  In  the  second  part  we  are  interested 
in  regularity  properties  of  6  in  r  and  need  estimates  that  are  uniform  for  V  G  ft,  r  G  B^.  Instead 
of  using  general  versions  of  the  implicit  function  theorem,  we  opt  for  a  direct  approach  that 
gives  precise  information  on  the  bounds.  Let  V  G  ft  be  fixed  and  consider  r,  s  G  B$  and  the 
corresponding  fixed  points  6{r)  and  6{s)  of  T.  Then  we  have 

(5.22)  Hr)  -  0(s)  =  [T(r,  9(r))  -  T( r, «(»))]  +  [Tf,  *(»))  -  T(s,  *(<))] . 


Using  (ii)  and  (iii)  in  Proposition  5.1,  we  obtain 


(5.23) 


Hence, 

(5.24) 


Mr)  -  0(«)||  <  ||T(r,  9(r))  -  T(r,  0(s))||  +  ||T(r,  0(s))  -  T(s,  0(«))|| 
<  a  ||0(r)  -  9(s) ||  +  C  6  |r  -  s| . 


and  (5.21)  follows  from  (5.6).  | 
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6.  The  map  connecting  the  wave  strengths  to  the  Riemann  data 

For  V  €  ft  the  states  V(£)  take  values  in  the  ball  B^U-)  =  {U  G  IRN  :  I U  -  U- \  <  fi}. 
Because  of  the  orthogonality  relations  (3.8)  and  the  continuity  properties  of  /,-({/)>  giyen 

T)>  0  we  can  choose  \i  such  that 

U(Ui)  ■  ri(U2)  >1-77,  UuUte  B„(U-) , 

(6.1) 

\U(Ui)  •  rj(U2)\  <  V ,  Ui,U2e  B^U-),  ifj. 

Also  for  states  in  B^{U-)  the  right  and  left  eigenvectors  are  bounded 

(6.2)  \ri(U)\  <  R,  \U(U)\<R ,  VeB^U.),  i  =  l,...,N, 

by  a  constant  R  depending  only  on  /i.  For  our  future  deliberations  we  place  an  additional  hypoth¬ 
esis,  which  complements  ( .1 1  —  Aj)  and  concerns  the  behavior  of  the  right  and  left  eigenvectors 
along  functions  in  ft:  Namely,  we  fix  T)  <  1/N  and  require  that 

(As) 

|/i(i7_)-ri(V(0)l  < 

for  V  €  ft  and  £  €  IR.  This  is  attained  by  restricting,  if  necessary,  the  size  of 
Consider  the  system  of  differential  equations 

N  N 

(6.3)  ea'k  +  [f- Xk(V(0)}  =  £  E  E  tVA*(F(0)  U^O)  •  ^(0)]  ««  “« , 

m=l  n=  1 

where  V  €  ft.  We  saw  in  the  previous  section  that  (6.3)  has  solutions  given  by  an  asymptotic 
expansion  in  a  parameter  r  €  of  the  form 

(6.4)  ak((  ;r)  =  Tk  <pk(()  +  9k({ ;  r) . 

The  expansion  is  valid  for  |r|  <  6  uniformly  for  V  €  ft,  and  <?*(•;  r )  satisfies  (5.20)  and  is  of  order 
0(|r|2)  as  |r|  — *■  0.  The  parameter  r  is  associated  with  the  data  at  ck ,  as  from  (5.1) 

(6.5)  ak(ck  ;  r)  =  Tk<pk(ck ) . 

It  is  instructive  to  visualize  |r|  =  |tj|  +  ...  +  |ryv|  as  measuring  the  wave  strength  of  the  Riemann 
problem  solution  associated  to  ak(£ ;  r)  (c./.  (3.11)). 

A  comparison  with  the  general  outline  in  Section  3  shows  that  while  the  solvability  of  (3.15) 
is  at  this  point  well  understood,  it  remains  to  select  r  so  that  (3.13)  is  satisfied.  The  issue  emerges 
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of  studying  the  connection  between  the  parameter  r  and  the  boundary  data  U±.  To  this  end  let 
V _  be  fixed  and  consider  the  map  S  that  carries  r  into  the  end-state  vector 

N  r  OO 

(6.6)  S(t)  =  U.  +  £  /  [r*  MO  +  r)]  rk(V(0)  d( . 

fc=i  J~°° 

For  t  £  Bs  =  {r  €  :  |r|  <  £}  the  map  S  is  well  defined,  and  depends  explicitly  on  V  and 

implicitly  on  e.  Our  objective  is  to  study  the  invertibility  of  5  and  to  show  that  the  inverse  map 
is  uniformly  bounded  in  V  and  e. 

Proposition  6.1.  Assume  that  ( A\  —  As)  are  satisfied  for  V  G  Cl.  There  exist  positive  constants 
r  and  6  such  that  : 

(i)  Given  U+  €  Br(U-)  there  exists  a  unique  solution  of  the  equation  S(r)  =  U+  with  r  £  Bs. 

(ii)  For  each  V  £  ft,  e  >  0  the  inverse  map  S-1  :  Br(U-)  —>  Bs  is  well  defined  and  satisfies 

(6.7)  \S-\U+)\<mU+-U.\, 

where  (3  is  a  constant  which  depends  on  fi,  but  is  independent  of  the  particular  V  €  and  £. 
Proof  of  Proposition  6.1.  Let  U-  be  fixed.  The  equation  5(r)  =  f/+  has  the  form 

N  r  OO  N  j.  oo 

(6.8)  U+-U.  =  Y/rk  Vk  rk(V(0)  <*C+£  /  MO  r)  rk(V(0)  d( . 

k=l  k=l  J~°° 

K  A(V)  is  the  matrix  whose  k-th.  column  is  given  by 

/oo 

-oo 

then  (6.8)  reduces  to 

N  ,oo 

(6.10)  U+-U.  =  A(V)  t  +  £  /  0k((;  T)  r*(F(0)  d( 

and  the  issue  becomes  to  study  the  solvability  of  (6.10)  in  r. 

First  we  show  that  Hypothesis  (-4s)  implies  that  A{V )  is  invertible. 

Lemma  6.2.  Assume  that  (A$)  holds  ( with  q  <  1/N).  The  matrix  A(V)  is  invertible  for  any 
V  £  Cl,  and  the  inverse  matrix  A~X{V)  is  uniformly  bounded 


(6.11) 


\A~\V)\<P,  V£Cl, 
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by  a  constant  (3  independent  of  e . 

Proof  of  Lemma  6.2.  Since  <pk  are  averaging  measures,  the  mean  value  theorem  implies 

,oo 

(6.12)  flfc(V)  =  /  <Pk  rk(V( 0)  =  rfe(Vfc*) 

for  some  V£  £  B^ ( U- ) .  Since  {r;(6r_ )}  are  linearly  independent,  by  choosing  fi  sufficiently  small 
it  is  guaranteed  that  the  vectors  ri(Fj*),  ... ,  r^r(V^)  are  linearly  independent  and  thus  A(V)  is 
invertible. 

We  now  show  (6.11)  and  in  the  process  provide  an  alternative  way  of  showing  that  A(V)  is 
nonsingular.  For  t,  y  £  IR^  consider  the  equation  A(V)r  =  y  and  write  it  in  the  form 


(6.13) 


roo 

YlTk  <PkTk(y(C))dC  =  y. 
k=\  J-co 

Taking  the  inner  product  of  (6.13)  with  li{U-)  and  rearranging  the  terms  we  obtain 

/oo  roo 

<Pi  [li(U-)  •  n(F(C))]  d(  =  h{u _)  ■y-J2T*  [W-)  •  d(  ■ 

-OO  i _ j.:  J  —oo 


k*i 


Then  (As),  (4.3)  and  (6.14)  yield 


(6.15)  N(1-t?)<  |/;(tf-)-S/|  +  »?2N- 

k*i 

Adding  the  resulting  equations  for  i  —  1,  ... ,  N  and  using  the  fact  that  rj  <  1/A,  we  obtain  the 
estimate 

1  N 

(6.16)  M  <  E I W-)  -v\<m  =  mv)r\. 

'  1=1 

The  first  implication  of  (6.16)  is  that  the  only  possible  solution  of  A(V)  r  =  0  is  the  trivial 
solution  t  =  0.  Therefore  ai(F),  ... ,  a;v( K)  are  linearly  independent  and  A(V)  is  invertible.  In 
addition,  (6.16)  implies  that 


(6.17) 


\A-1(V)y\</3\y\,  ye  IR", 


which  proves  (6.11).  | 

Next,  we  formulate  solving  the  equation  S(t)  =  U+  as  a  fixed  point  problem.  Let  Br(U-) 
be  the  ball  centered  at  U-  of  radius  r  and  consider  the  map  P  that  takes  U+  €  Br(U-),  V  £  Q, 
t  £  B$  into  the  vector 

N  foo 

(6.18)  P(U+,V,t)  =  A-\V){U+-U-)-A-\V)Ys  /  «*«;  r)r*(V(C)K- 

k=  iJ-°° 
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Since  A(V )  is  invertible,  solutions  of  (6.10)  are  equivalent  to  fixed  points  of  the  map  P(U+ ,  V,  •). 

Lemma  6.3.  There  exist  positive  constants  6  and  r  such  that  P  :  Br(U-)  x  ft  x  Bg  — ►  Bg  and  has 
the  property  :  There  exists  a  constant  a  with  0  <  a  <  1  such  that 


(6.19)  \P(U+,V,T)-P(U+,V,s)\<a\r-s\,  r,  s  €  Bs , 

for  any  U+  €  Br(U _),  V  6  fi;  that  is  P(U+,V,-)  is  a  uniform  contraction  on  Bg. 

Proof  of  Lemma  6.3.  Let  U+  €  Br(U _),  V  €  Q  and  r  €  Bg.  Using  (6.18),  (6.11),  (6.2)  and 

(5.20) ,  we  obtain 

\P(U+,V,t)\<\A-\V)\  (\U+-U. |  +  S/I  I4(C;  r)|MU(0)|dc) 

(6.20)  /  JL  r°°  \ 

<(3{r  +  RA\T\2NY,J~  <pjd^ 

<  /3(r  +  RAN2  S2) 

The  first  part  of  the  lemma  is  true,  provided  that  r  and  6  satisfy 

(6.21)  pr  +  PRAN2  62  <8. 


Let  now  r,  s  6  Bg  and  observe  that 

N  too 

(6.22)  P(U+,V,t)~  P(U+,V,s)=  -A-^V)^  /  [9k«;  r)  -  9k((;  *)]  rk(V(0)d(. 

k= iJ-°° 

On  account  of  (6.11),  (6.2)  and  (5.21),  (6.22)  gives 

N  foo 

I  P(U+,  V,T)-P(U+,V,s)  |  <  /?  £  /  l^(C;  r)  -  0k(C ;  s)\  |rfc(V(C))|  d( 

k= lJ-°° 

N  .oo 

<0RNC8\t-s\  J2  /  <Pjd( 

1  j  —  oo 


(6.23) 


j- 1 


Therefore,  if 


(6.24) 


<  p  R  N2  C  8  |r  -  s\ . 


a  —  (i  RN2  C  8  <  1 


then  P(U+,V,  •) :  Bg  — ►  Bg  is  a  uniform  contraction. 
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Note  that  if  6  <  Amin {(/LRJV2C)-\  (PRN2A)~1}  and  r  <  ±6  then  both  (6.21)  and  (6.24) 
are  simultaneously  satisfied,  and  the  proof  of  the  lemma  is  complete.  | 

We  return  to  the  proof  of  Proposition  6.1.  Lemma  6.3  implies  that  given  U+  €  Br{JJ-)  there 
exists  a  unique  fixed  point  of  P(U+,V,‘)  in  the  ball  B$  and  thus  a  unique  solution  of  S(r)  =  {/+. 
Hence,  S~*  is  well  defined.  Let  {7+  and  r  =  5-1(i 7+)  be  two  corresponding  points  connected 
through  (6.10).  Using  (5.20),  (6.2)  and  (4.3),  we  obtain 

N  roo 

\A(V)t\<\U+-U-\  +  Y,  /  \0k«;r)\\rk(V(O)\d( 

(6.25)  ,  A  f°° 

<\U+-U-\  +  RAN\t\2J2  tpjdt 

j= 1  J-°° 

=  \U+-U.\  +  RAN2\t\2 

Using  Lemma  6.2,  in  conjunction  with  (6.21)  and  the  choice  of  6 ,  we  deduce  from  (6.25) 

(6.26)  |r|  </3\U+-U-\  +  PRAN26\t\<P\U+-U-\  +  ^\t\, 
which  implies  (6.7)  and  completes  the  proof  of  the  proposition.  | 
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7.  Proof  of  Theorem  3.1 


This  is  the  concluding  section  of  the  derivation  of  a-priori  estimates  for  (Vc).  The  analysis  of 
Sections  3  to  6  is  combined  in  order  to  prove  the  main  theorem. 

Let  U-  be  fixed  and  define  ft  by  (4.1).  ft  is  a  convex,  bounded  subset  of  the  Banach  space 
C°(— oo,  oo)  of  continuous,  bounded  functions.  Fix  e  >  0  and  consider  the  map  T  carrying  V  €  ft 
to  the  continuous  function  W  defined  by  the  procedure  : 

(a)  Let  (pk  be  as  in  (4.3).  We  obtain  the  solution  0*(-;r)  of  (5.1),  for  r  €  IR^  small,  and  define 

<**(•  5 T)  =  Tk<fik  +  ;  t).  The  resulting  afc  form  a  solution  of  the  system  of  equations  (6.3). 

(b)  Let  S  be  the  map  defined  in  (6.6).  Let  t  be  the  solution  of  the  equation  S(t)  =  U+,  that  is 
t  =  S-\U+). 

(c)  W  is  then  defined  by  setting 

ri  N 

(7.1)  w(0  =  u.+  X  [**  MO  +  **(C;  t)}  n(v(0)  d( . 

J-°°  k= i 

Concerning  the  feasibility  of  the  constructions  we  remark  :  The  parameter  p  in  the  definition 
of  ft  is  fixed  so  that  Hypotheses  (A\  —  As)  are  satisfied  for  V  6  ft-  We  also  fix  the  parameters  A 
and  So  as  in  Proposition  5.1  and  let  6  <  So.  Then  Corollary  5.2  states  that  for  r  6  Bg  the  problem 

(5.1)  has  a  unique  solution  satisfying  the  estimate 

N 

(7.2)  \Ok(-;r)\  <  A\t\2  t  6  Bg . 

i= i 

According  to  Proposition  6.1,  for  r  and  6  sufficiently  small  the  map  S  :  Bs  — ►  Br(U-)  is  invertible, 
S(t)  =  U+  is  uniquely  solvable  in  Bg,  and  the  inverse  t  =  S'- 1  (?7+)  satisfies  for  some  fixed  (5 
(independent  of  V  and  £•)  the  estimate 

(7-3)  |*|  =  \S~\U+)\  <20\U+-U-\,  U+<=  Br(U-) . 


As  a  result  W{— oo)  =  U- 
(7.4) 


and  W(+oo)  =  S(t )  =  f/+.  From  (7.1)  we  obtain 

^-  =  X  [**?*  + **(•;*)] 


which,  in  conjunction  with  (7.2)  and  (6.2),  yields 


(7.5) 


l^-l  ^  RY,  [lt*IP*  +  4*l3£v,i]  |r*(V)| 

N 

<  -R|*|  (1  +  AN  |*|)  X  Vj  • 
i- 1 
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In  turn,  (7.1),  (7.3)  and  (7.5)  imply 


(7.6) 


|w(£)-£/_|<|  I  £[frn(C) +  **«;<)]  rk(v(Q)\dc 

\J-oo  k=l 

<  2 pNR  (1  +  2/3AN  \U+  -  tf_|)  \U+  -  V-\ . 


It  follows  that,  if  U+  €  Br(U _)  and  r  is  restricted  by 


(7.7) 


2PNRr(l  +  2pANr)  </z, 


the  function  W  defined  in  steps  (a  -  c)  satisfies 


(7.8) 


\W(i)-U-  \</i,  {€1R. 


In  the  sequel  we  fix  r  and  S  to  simultaneously  satisfy  (7.7),  (6.21),  (6.24)  and  (5.16).  All  the 
stated  constructions  and  estimations  are  then  feasible,  and  the  map  T  :ft  ^  ft  is  well  defined.  In 
addition,  (7.5),  (7.3)  and  Lemma  4.1  dictate  there  is  a  constant  C  such  that 


(7.9) 


\W(0~U-\<\U+-U-\-  f*  e-£(c-Al-)2dC  7  for  £  <  Ai_  , 

£  J -oo 

C1  t°° 

\W(Z)  -  U+ 1  <  \U+  -  U-l  -  /  e-*(c-Aw+)  d(  ,  for  £  >  \N+  . 


Our  next  task  is  to  apply  the  Schauder  fixed  point  theorem  to  the  map  T. 

(i)  T(Q)  is  precompact  in  C°(— oo,  oo). 

Consider  a  sequence  {Vn}  C  ft  and  let  Wn  =  T(Vn).  Estimates  (7.5),  (7.3),  (7.8)  and  (4.7) 
imply  that  { Wn }  is  uniformly  bounded  and  uniformly  equicontinuous  on  the  reals.  It  follows  from 
the  Ascoli-Arzela  theorem  and  a  diagonal  argument  that  there  is  a  subsequence  {Wnj }  and  a 
continuous  function  W  such  that  Wn*  — ►  W  uniformly  on  compact  subsets  of  IR.  But  then  the 
decay  estimates  (7.9)  imply  that  the  convergence  is  in  fact  uniform,  and  thus  T(ft)  is  precompact 
in  C°(— oo,  oo). 

(ii)  T  :  ft  -*■  ft  is  continuous. 

Let  {Vn}  C  ft  be  a  convergent  sequence  in  C°(— 00,00),  with  Vn  — +  V°,  and  set  Wn  =  T(Vn), 
W°  =  T(V°).  We  proceed  to  show  T(Vn)  T(V°).  RecaU  that  e  is  held  fixed,  and  that  Wn 
and  W°  are  defined  in  terms  of  the  intermediate  quantities  <pk,  9k(-  ;r),  ak(-',T),  Sn,  tn  and  <p ”, 
<?2(* ;  t)  a°k( • ;  r),  S°,  t°  in  steps  (a  -  c)  for  V  =  Vn  and  V  =  V°,  respectively. 

First,  we  show  <pk  — »•  <p°k  in  C°(- 00, 00).  One  first  uses  (4.7),  (4.2)  and  (4.8)  to  show  that  {<^£} 
is  a  uniformly  bounded  and  equicontinuous  sequence  of  functions  that  satisfies  the  decay  estimates 
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(4.8)  as  |£|  — ►  oo.  An  argument  as  in  (i)  implies  there  exists  a  subsequence  {<pk’ }  and  a  function 
<pf  such  that  <£>"'  — *■  uniformly  in  1R.  Parsing  to  the  limit  in  (4.3)  along  the  subsequence  rij 


(7.10) 


w  = 


fO°  e~‘  ft  (»))  d»  dc  roo  e-j  f'  s-Xt(V0(3)) 

J— OO  ’  J  —  oo 


ds 


= 


dc 


we  deduce  (p =  <p°k.  The  sequence  {y>£}  has  limit  points,  and  any  limit  point  is  equal  to  (p°k. 
Hence,  the  whole  sequence  {<£>£}  converges  to  <p°k. 

Second,  we  show  that  for  r  fixed  6k( • ; r)  —*■  6°k(- ;  r)  in  C°(— oo,  oo).  This  follows  by  a  similar 
in  spirit  argument  that  we  only  sketch:  Using  (5.20),  (5.1),  (4.7)  and  (4.8),  we  show  that  {<?£} 
possesses  a  subsequence  { 9k 3 }  and  a  limit  point  0£°  so  that  6k’  — *•  0£°  uniformly  in  1R.  Passing  to 
the  limit  in  (5.4)  along  the  subsequence  rij  and  using  the  convergence  of  Vn  and  <pk,  we  obtain 


(7.11)  6?(0  = 


=  6  « 9k  I  e*9k 
Jck 


N 

m,n= 1 


.(V°(0)  (r.  A(0  +  « 


(0)  (r.vt(0  +  C(0)<fC. 


Since  the  limiting  6 £°  inherits  the  estimate  (5.20),  the  uniqueness  part  of  Corollary  5.2  implies  any 
limit  point  of  {#£}  is  of  the  form  0£°(- ;  r)  =  9°k(  - ;  r).  Consequently  0k(- ;  r)  — ►  6°k(  - ;  r). 

The  third  step  is  to  show  that  tn  — *•  t°  in  1RN.  Let  5"  and  S°  be  the  maps  associated  with  Vn 
and  V°  respectively  and  define  tn  and  <°  satisfying  Sn(tn)  =  S0^0)  =  {/+.  Since  {<"}  is  bounded, 
there  is  a  subsequence  {tn> }  and  a  vector  t°°  such  that  tn>  — ►  t°°.  We  use  (5.20),  (5.21)  to  pass  to 
the  limit  in  Sn>(tn> )  =  U+  and  to  obtain  S0(t°°)  =  U+.  Because  of  the  unique  invertibility  of  the 
map  S°  it  is  t°  =  t°°  and  thus  the  sequence  {fn}  converges  to  t°. 

The  precompactness  of  T  implies  the  sequence  {Wn}  has  a  subsequence  {Wn> }  and  a  limit, 
function  W°°  such  that  Wn>  — *•  W00  in  C°(- oo,oo).  Using  the  established  convergences  and 
(5.21),  we  pass  to  the  limit  in  (7.1)  along  nj  and  obtain 


(7.12) 


W°°(0  =  U-  + 


[<*  vS(C)  +  0”(C;  t °)1  r*(v°(0)  dc  =  T(V°)(  0 . 


Therefore  any  limit  point  of  {Wn}  is  equal  to  T(V°)  and  thus  T(Vn)  T(V° )  in  C°(-oo,oo). 

Hence,  T  is  continuous. 

The  Schauder  fixed  point  theorem  implies  there  exists  a  fixed  point  Ue  of  the  map  T  in  fL  By 
construction  JJe  satisfies 


(7.13) 


uc(0  =  U.+ 


N 

2Jake((;Te)n(ve(())  d(, 


k- 1 
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where 


(7.14) 


akc(£'i  Tc)  —  Tk,c  V’fceCO  +  ^fce(£)  Te) 


solves  (3.15).  The  functions  <pke,  0ke  and  ake  depend  implicitly  on  e,  and  the  quantities  r£  satisfy 

N  fco 
k=  1  •,-°° 

As  a  result  U£(±oo)  =  U±  and 


(7.15) 


N  roo 

5£(t£)  =  U-  +  Y,  /  te)  rk(UE( 0)  d(  =  £/+ 

I _ 1  J -CO 


N 


(7.16) 


U'c(0  =  ^2  ake({;  Te)rk(Ue(0) , 

Ar=l 

«*.«;  r.)  =  ))•«({). 

Using  (7.16)  and  (3.5  —  3.8),  we  can  rewrite  (3.15)  in  the  form 


(7.17) 


lk(Uc)  •  [  -  f  +  VF(t/£)  ]  U'e  =  lk(U£ )  •  U'J 


which  implies  that  f/£  is  a  solution  of  (Pc). 

Consider  a  family  {U£}£>o  of  such  solutions  to  (P£).  By  construction  Ue  are  of  uniformly 
bounded  (and  small)  oscillation  ( C0 )  and  satisfy  the  representation  formula  (3.33).  Relations  (7.2) 
and  (7.3)  imply  there  exist  constants  C  independent  of  £  so  that  |r£|  <  C\U+  —  U- 1, 


(7.18) 


N 


|flfce(£;  r£)|  <  \TktC\(pke  +  C|r£|2 


3= 1 


N 


<  C \U+  -  U- \{<fike  +  | u+  -  U-|  £>;,) 

3= 1 


and  as  a  result 
(7.19) 


AT 


|t%)l  <  K  £>* 

i= i 


with  K  a  constant  of  order  0(|{7+  —  U- 1)  and  independent  of  e.  As  {ipj£}  are  uniformly  bounded 
in  X1(IR),  it  follows  that  {U'e}  are  uniformly  bounded  in  Ll  (IR)  and  {U£}  is  of  uniformly  bounded 
variation.  The  total  variation  of  the  family  is  controlled  by  \U+  —  U- 1  and  is  thus  small.  The  proof 
of  Theorem  3.1  is  complete. 
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8.  Solution  of  the  Riemann  problem 


Our  next  objective  is  to  construct  solutions  of  the  Riemann  problem  (V)  by  taking  e  —*  0 
limits  of  solutions  of  ( Vc )  and  to  identify  the  structure  of  the  emerging  solutions.  The  analysis  is 
patterned  within  the  framework  developed  in  the  previous  sections.  Nevertheless,  it  is  instructive 
to  single  out  the  set  of  hypotheses  used  in  performing  the  £  — ►  0  limit  and  to  provide  an  independent 
presentation.  Let  {t/£}E>o  be  a  family  of  solutions  to  ( Ve )  that  connect  U-  to  U+  and  enjoy  the 
properties  : 

Ue  satisfy  the  uniform  bounds  (C0),  ( S )  for  e  >  0  , 

( Aa )  ^k(Uc)  satisfy  the  uniform  bounds  (3.9),  (3.10)  for  e  >  0  , 

U'e  satisfies  (7.19)  where  (fkc  is  given  by  (3.25) . 

Solutions  satisfying  (/t^)  were  constructed  in  Theorem  3.1,  and  the  resulting  families  are  of  small 
oscillation  and  variation.  The  results  of  this  section  remain  valid  for  families  of  large  oscillation 
and  variation,  provided  the  global  separation  of  the  eigenvalues  and,  most  important,  estimate 
(7.19)  hold.  Helly’s  selection  principle  implies  there  exists  a  subsequence  of  the  original  family, 
denoted  again  by  {Uc}  with  £  — ►  0,  and  a  function  U  of  bounded  variation  such  that 

(8.1)  t/e(f) -*■  {/(£)  pointwise  on  (-00,00). 

Since  U  is  of  bounded  variation  its  domain  can  be  decomposed  into  two  disjoint  subsets  C  and  S  : 
C  consists  of  the  points  of  continuity  of  U  and  S  of  the  points  of  jump  discontinuity.  S  is  at  most 
countable,  and  the  right  and  left  limits  of  U  exist  at  any  £  €  S  and  are  denoted  U(£±). 

We  proceed  to  show  U  satisfies  ( V ).  In  the  sequel  C  denotes  a  generic  constant  that  can  be 
estimated  in  terms  of  the  bounds  in  (/ls)  and  the  Riemann  data  and  which  is  independent  of  e. 

Theorem  8.1.  Suppose  that  (1.1)  is  strictly  hyperbolic  and  let  {(/£)}5>o  be  a  family  of  solutions  of 
(Ve)  corresponding  to  data  U±  and  satisfying  (j4s).  There  exists  a  subsequence  {UCn}  with  en  -*  0 
and  a  function  of  bounded  variation  U  such  that  U£n  — ►  U  pointwise  on  the  reals.  U  satisfies 

(8.2)  -tU'  +  F(U)'  =  0 

in  the  sense  of  measures ,  the  Rankine-Hugoniot  conditions  hold  at  any  point  £  €  S 


(8.3) 


-e  Mt+)  -  U(t-)]  +  [F(U(t+))  -  F(U(f-)]  =  0  , 


and  there  exist  constant  vectors  Uq,  ,  Un  G  with  Uo  =  U-,  —  U+  such  that 

'  Uq  =  U-  -oo<e<A1_, 

(8.4)  U(0  =lUk  Xk+  <  £  <  \k+i)-  ,  k  =  1,  ... 

k  Un  =  U+  \n+  <  £  <  +00  • 

Proof.  Let  {U£}  be  a  convergent  subsequence  as  in  (8.1),  satisfying  the  uniform  bounds  (C0),  {S), 
and  ip  €  [(^(IR)]^  be  a  test  function  with  compact  support.  Then  ( V£ )  gives 

(8.5)  [  U£  •  (£V>)'  -  F(Uc)  -1>'dt  =  e  f  U£-  d( . 

J IR  ^IR 

Passing  to  the  limit  £-+0we  deduce 

(8.6)  /  U-{W)'-F(U)-i>'dt  =  0 

J  IR 

that  U  satisfies  (8.2)  in  the  sense  of  distributions.  Since  U  is  of  bounded  variation,  it  also  satisfies 
(8.2)  in  the  sense  of  measures. 

Let  £  =  [Ai_,  Ai+]  U  ...  U  [A^_,  Ajv+]  stand  for  the  range  of  variation  of  the  wave  speeds 
A k(Uc).  Then  (4.8)  and  (7.19)  imply 

C  1 

<Pkc  <  y  exp{-— d(£,Afe)2}  £  €  (-00,00)  -  [A*_,  Afc+] , 

(8.7)  n  r  1 

|CA'|<if^^£<^exp{--da,£)2}  £  €  (—00, 00)  —  C , 

j- 1 

where  d(£,  \k)  and  d(£,£)  are  the  distances  between  the  point  £  and  the  sets  [Afc_,  A*+]  and  C 
respectively.  Therefore  the  limiting  function  U  stays  constant  on  each  connected  component  of 
(—00,00)  -  C  and  (8.4)  follows.  In  addition  t/£(±oo)  =  U±  implies  Uo  =  U-  and  t//v  =  U+. 

The  Rankine-Hugoniot  conditions  (8.3)  are  a  consequence  of  the  fact  that  U  of  bounded 
variation  solves  (8.2).  We  outline  a  different  proof,  in  the  spirit  of  self-similar  viscous  limits. 
Integrating  the  equation  (V£)  on  an  interval  (a,  b)  we  obtain  the  weak  form 

rb 

(8.8)  [-b  U£(b)  +  F(U£{b ))]  -  [-a  U£(a )  +  F(tf£(a))]  +  /  U£(()  d(  =  eU'£(b)  -  eU'e{a) . 

J  a 

For  f  €  S  and  6  >  0,  we  evaluate  (8.8)  between  the  points  0  and  r,  with  r  <  Ai_,  and  integrate 
the  resulting  equation  in  9  over  [f ,  £  +  6]  to  arrive  at  the  identity 


-0U€{0)  +  F(Ue(0))  d0  +  j  J  Ue(()d(d0 

=  £  J  U'e{0)d0-e6U'E{T)  +  6{-rUe{T)  +  F{U£{T))\. 


(8.9) 
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Using  the  consequence  of  (7.19)  and  (3.25) 


(8.10) 


/oo 

\U'c\d(  <  KN  , 

-OO 


in  conjunction  with  (A„),  (8.1)  and  (8.7),  we  take  first  e  -+  0  in  (8.9)  and  then  divide  the  resulting 
equation  by  6  and  take  6  -*•  0+  to  obtain 


(8.11) 


-(U((+)  +  F(U((+))  +  £  V(()d(  =  -r  V(t)  +  F(V(t))  . 


In  a  similar  manner,  given  any  6  and  r  <  Ai_,  we  establish 

(8.12)  -eu(0-)  +  F(U(0-))  +  J$  U(0d(  =  -T  U(t)  +  F(U(r))  . 

Then  (8.3)  follows  from  (8.11)  and  (8.12)  for  £  =  6.  | 


With  U(£)  as  above,  define 

(8.13)  V(x,t)  =  U(y),  (x,t)  e  (-00,00)  x  (0,  oo ) 

Clearly  lim^o  V(x,t)  =  {/_  for  x  <  0,  U+  for  x  >  0.  Furthermore,  a  solution  V  of  the  form  (8.13) 
is  a  weak  solution  of  (1.1)  on  (—00,00)  X  (0, 00)  if  and  only  if  U  is  a  weak  solution  of  (8.2)  on 
(-00,00).  The  equivalence  follows  from  an  argument  due  to  Dafermos  [D3].  Let  x(z,f)  be  a  C°° 
IR^- valued  function  with  compact  support  in  (-00,00)  x  (0, 00)  and  define 

f°° 

(8.14)  V>(£)  =  /  x(£M)*« 

The  resulting  function  ip  €  [C^°(— 00, 00)]^.  Conversely,  any  test  function  ip  may  be  represented 
in  the  form  (8.14)  by  choosing  x  =  V’(x/0a(0>  °(0  €  C'^°(0,oo)  a  fixed  function  such  that 

/0°°  a(t)  dt  =  1.  For  solutions  of  the  type  (8.13)  the  weak  form  of  (1.1)  may  be  written  as 

[CO  [OO 

/  /  V(x,t)-Xt(x,t)  + F(V(x,t))-Xx(x,t)dxdt 

Jo  J-00 

/OO  [OO  [OO 

Xt(tit,t)tdt)  +  F(U(  £))•(/ 

-00  Jo  Jo 

/oo 

-00 

and  the  equivalence  follows  from  the  chain  of  identities.  Theorem  8.1  in  conjunction  with  Theorem 
3.1  lead  to  an  existence  theorem  for  the  Riemann  problem. 
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Theorem  8,2.  Assume  that  (1.1)  is  strictly  hyperbolic .  Given  any  data  U+  with  \ U+  —  U-\ 
sufficiently  small ,  there  exists  a  function  of  bounded  variation  U (£)  defined  on  (—00,00)  such  that 
C7(  y  )  is  a  weak  solution  of  the  Riemann  problem  for  (1.1). 

Next  we  investigate  the  structure  of  the  emerging  solution  U .  It  is  instructive  to  use  the 
correspondence  between  functions  of  bounded  variation  and  finite  signed  Borel  measures  on  IR 
(Folland  [F,  Sec  3.5,  Sec  7.3]).  Let  p  be  the  (vector  valued)  measure  generated  by  the  right 
continuous  function  of  (normalized)  bounded  variation  ({/(£+)  ”  U-).  Consider  now  the  functions 

(8.16)  <M 0=  f  <Pk*(QdC. 

J —00 

In  view  of  (3.25)  the  family  {$jte}  consists  of  increasing  uniformly  bounded  functions.  Therefore 
converge  along  a  subsequence  to  an  increasing  function  point  wise  on  the  reals.  The  measures 

generated  by  $*(£+)  are  denoted  by  4 ;  they  are  positive  measures  with  total  mass  one. 
Introduce  the  measures  associated  with  the  functions  U'e  and  (fke  defined  by 


(8.17) 


</z£,V>>  =  [  U'£(0-mdt, 

J  IR 

<  4e>X  >  =  f  (Pke(0x(Qd£, 


where  xp  G  [C^IR.)]^,  X  €  Cc(IR)  are  continuous  functions  with  compact  support.  Then  (3.25), 
(7.19),  (8.10)  and  Helly’s  convergence  theorem  imply 


(8.18) 


/  u'e-«>d(^  [  xp-du  =<  n,xp> ,  forv>e  [cc(nt)]N, 
J  IR  J  IR 

I  VkcXdt.^  I  xd$k  =<  <t>k,x  >,  for  X  €  Cc(IR) . 

JlR  J  IR 


In  the  language  of  functional  analysis  /ze  — 1 ■  n  and  4e  — ' ■  4>k  weak-*  in  measures. 

Using  (8.18)  we  can  express  <  /z,  xp  >=  —  U  •  xp'd£  for  test  functions  xp  6  [^(IR)]^.  Note 
that  f  £  supp/z  if  and  only  if  there  is  an  open  interval  I  B  (  such  that  <  //,  xp  >=  -  =  o 

for  xp  6  [C(!(/)]Ar.  This  is  in  turn  equivalent  to  the  function  U  being  a.e.  equal  to  a  constant  vector 
on  I.  Consequently  supp  fi  coincides  with  the  region  in  the  ^-domain  where  U  is  not  a  constant 
state.  From  (8.7)  it  follows  that  /z  is  absolutely  continuous  with  respect  to  $k  and  that 


(8.19) 


supp  4  C  [A*_,  Afc+] 

N  N 

supp  /z  C  (J  supp  4  c  C  =  (J[Ak_,  Ak+]. 

k=  1  k=l 
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The  following  proposition  states  an  important  property  of  <j>k ,  that  incorporates  admissibility 
restrictions  induced  by  the  self-similar  viscosity.  In  preparation,  recall  that 


(8.20) 

where 

(8.21) 


<fke 


e~i3k‘ 


SZ>e~‘9k‘dC 


0fce(£)  —  f  S-Xk(U£(s))ds, 

J  Pkt 


and  assume  (by  restricting  to  a  further  subsequence)  that  pkc  — >  pk  as  e  — *  0.  Using  (8.1),  (C0) 
and  the  Ascoli-Arzela  theorem  we  deduce  that 

(8.22)  gke(0  =  f  s  -  A k{Uc(s))  ds  — >  f  s-  A k(U(s))  ds  =:  gk( £) 

Jpki  •'Pk 

uniformly  on  compact  subsets  of  (—00,00).  We  show  that  points  in  the  support  of  <f>k  are  global 
minima  for  the  function  gk. 

Proposition  8.3.  //£  G  supp<^  then  gk(0  >  9k{0  for  C  £  (—00,00). 

Proof.  The  proof  has  two  steps.  First,  fix  any  £  €  IR,  and  a  >  0  and  consider  the  set 


(8.23)  A  =  (C  €  IR  :  fffc(C)  ~  9k(0  <  -a  <  0} 

Since  gk  is  continuous  either  A  is  empty  or  it  has  positive  Lebesgue  measure  m(A).  We  will  prove 
Claim :  If  m(A)  >  0  there  exists  an  open  interval  /  9  £  such  that  <  (f>k,X  >=  0  for  any  X  £  Cc{I). 

As  a  result,  if  m(A)  >  0  then  £  £  supp  <fo. 

To  establish  the  claim,  observe  first  that 

(8.24)  flf*(C)  "  9k(0  >  \(C2  ~  e2)  -  max{|A*_|,  |Afc+|)  |C  -  £| 

implies  jffc(C)  — ►  00  as  |£|  — >  00  and  A  is  contained  in  some  compact  interval  [a,  b}.  Fix  >  0  such 
that  for  9  G  (£  -  tf,£  +  6)  we  have  |srfc(£)  -  </*(0)|  <  f .  By  virtue  of  (8.22),  there  is  e0  >  0  such 
that  if  £  <  £0  then 

(8.25)  \9ke(6)  ~  9k(6)\  <  ^  for  0  G  A  U  (£  -  S,  £  +  6). 

From  (8.23)  and  (8.25)  we  deduce,  if  0  G  (£  -  £,£  +  £),  £  <  £0  and  C  £  >4  then 

9kc(C)  ~  9kc(0)  <  gk(0  ~  9k(0  +  1 9k(0  ~  9k(0) I 

^  ^  +  \9ke(0)  ~  0fc(0)l  +  IflMC)  -  5fc(C)l  <  • 
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In  turn  (8.20)  and  (8.21)  yield  for  9  6  /  :=  (£  -  +  £) 


(8.27) 


1 

o  <  <pke{0)  <  j^exp{-i(fffce(0  -9ke(0))}d(  -  m{A) 


Let  x  G  Cc(I).  Then  (8.18)  and  (8.27)  imply 


(8.28) 


<<t>ke,X>=  f  <Pkc{O)x(0)M-+  °>  as  e-^0. 


Hence,  <  <f>k,X  >=  0  for  x  G  Cc(I),  and  the  proof  of  the  claim  is  complete. 

Suppose  next  that  £  €  supp<£k.  Then  A  is  empty  for  any  a  >  0  and,  hence,  gk(0  >  ffk(0  for 
any  £  £  (— oo,  oo).  | 

The  minimization  properties  for  the  gk  yield  information  on  the  structure  of  U .  In  particular, 
a  weak  form  of  the  Lax  shock  conditions  is  induced  at  points  of  discontinuity. 

Proposition  8.4.  Let  £,  £'  £  supp  p  fl  [Ak_,  Ak+]  with  £  <  £'. 

(a)  If  £  £  C  then 


(8.29) 


£  =  h(u(0)  • 


(b)  //£  €  5  then  U  satisfies  at  £  the  jump  conditions  (8.3)  and  the  inequalities 


(8.30) 


A k(U(t+))  <  £  <  A k(U(t-))  ■ 


(c)  //£,  £'  G  supp  p  n  [Ak_,  Ak+]  then  A *(£/(£+))  =  £,  A *(£/(£'-))  =  £'.  Moreover,  for  any  point 

*€(£,£') 


(8.31) 


9  =  X  k(U(0))  if  6eC, 
Xk(U(6+))  =  0  =  Xk(U(0-))  if  oes. 


Proof.  The  function  gk  in  (8.22)  is  continuous  and  has  the  behavior  gk(0  00  as  l£|  -+  00 .  Since 
U  is  of  bounded  variation,  the  limits 


(8.32) 


lim  gfc(C)  -  gfc(0 
c-*±  <  -  £ 


=  lim  — [  s  -  Xk(U(s))ds  =  £  -  Xk(U(£±)) 
C->£±  (,  -  $  Je 


exist  and  imply  the  derivative  exists  and  is  continuous  for  (  f  C,  while  only  the  right  and 
left  derivatives  exist  for  £  G  S.  Fix  a  point  £  £  supp  p  D  [Ak_,  Ak+].  It  follows  from  (8.19)  and 
Proposition  8.3  that  £  £  supp  <^>k  and  that  gk(()  >  gk(0  f°r  C  £  IR-  In  turn,  (8.32)  yields 


(8.33) 


£-Afc(tf(£+))>0  ,  £-Afc(C7(£-))<0, 
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which  leads  to  (8.29)  for  £  £  C  and  to  (8.30)  for  £  £  S. 

It  remains  to  show  (c).  Let  £,  £'  €  supp  fi  fl  supp  fa  with  £  <  £'.  Then  £,  £'  are  both  global 
minima  for  gk  with  <?*(£)  =  9k(t')-  We  claim: 

(8.34)  9k(0)  =  gk(O  for  any  6  £  (£,£') . 

If  (8.34)  is  violated  at  some  point,  there  exist  a,  b  with  £  <  a  <  b  <  £'  such  that 

(8.35)  gk(a)  =  gk(b)  =  gk(£)  ,  gk(0)  >  gk(£)  for  a<9  <b. 

At  the  points  a,  b  we  have 

A  k(U(a+))  <a<  Xk(U(a-)) 

(8.36) 

Xk(U(b+))  <  b  <  A k(U(b-)) . 

On  the  other  hand  at  any  6  £  ( a,b )  the  set  A  =  {C  G  IR  :  gk(()  —  gk{9)  <  — »}  is  nonempty  for 
some  a  >  0.  Proposition  8.3  and  (8.19)  then  imply  0  ^  supp  fa  and  the  function  U (£)  remains 
constant  on  the  interval  (a,  b).  Hence  A^(f7 (a+))  =  \k(U(b-))  and  the  inequalities  (8.36)  yield 
b  <  a.  This  contradicts  a  <  b  and  (8.35)  follows.  | 

In  summary,  the  region  where  U  is  nonconstant  consists  of  (at  most)  N  disjoint  closed  intervals 
I\k  =  [ak,  &*],  k  =  1,  ... ,  N.  Each  I\k  is  associated  with  one  characteristic  speed  \k(U)  and  could 
be  empty  or  consist  of  just  a  single  point.  The  function  U  takes  constant  values  on  the  complement 
of  UfcLi  and  has  the  properties  listed  in  Proposition  8.4  at  points  of  I\k .  The  emerging  solution 
consists  of  k- wave  fans  separated  by  constant  states.  Next  we  use  the  weak  form  of  (8.2) 

(8.37)  -£  U(£+)  +  0  U(0- )  +  F(t/((£+))  -  F(U(0-))  +  f  U(s )  ds  =  0  £,  0  £  IR . 

Je 

in  conjunction  with  relations  (8.29  —  8.31)  to  obtain  a  fuller  description  of  the  behavior  of  U  on 
the  wave  fans. 

Proposition  8.5  Suppose  that  I\k  =  [ak,bk]  is  a  full  interval,  ak  <  bk. 

(i)  For  each  £  6  [a*A)  such  that  VA *([/(£+))  ■  rfc(t/(£+))  ±  0 

(M8)  » -*>«  m+))  -  nm+)) 

(ii)  For  each  £  6  (ak,bk]  such  that  VA *(£/(£-))  •  rk(U({-))  ?  0 

hJomh<o  h  +  h+ )  ”  =  VA*(If(£-))-rfc(17(£-)) 


(8.39) 


Proof.  We  show  (i).  Fix  £  €  [a*,6fc)  and  let  h  >  0  such  that  £  +  h  €  I\k-  The  weak  form  (8.37) 
taken  between  the  points  £+  and  £  +  h—  gives 

[  -  (/  +  VF((/(?+))  ]  (P({  +  />-)  -  V((+)) 

(8.40)  =  -  [  W((  +  ft-))  -  F(U((+))  -  VF(U  (,(+))  (V((  +  ft-)  -  Plf +))  ] 

-  j  [  J/(4j)  -  £/(£+)  ]  *  +  (tf(£  +  h-)  -  £/(£+)) 

The  increment  (£/(£  +  h— )  —  !/(£+))  is  expanded  in  the  basis  of  right  eigenvectors 


(8.41) 


u(h)  :=  t/(£  +  h-)  -  t/(£+)  =  5>;(/i)  r,(tl(£+)) . 


Note  that  for  a  function  U  of  bounded  variation  u(h )  ->0  a s  h  -+  0+,  and  that  by  (3.8) 


(8.42) 


Ui(h)  =  U(U(t+))-u{h). 


Taking  the  inner  product  of  (8.40)  with  /;(#(£+))  and  using  (3.6),  (8.42)  and  the  Taylor  expansion, 
we  obtain 

(8.43)  [ -£  +  A, •({/(£+)) ]  u>i(h)  =  0(|u;(h)|2)  +0{j‘  \u(s)\ds)  +  0(hMMI) , 

On  account  of  Proposition  8.4  and  the  strict  hyperbolicity  of  (1.1),  the  coefficient  [— £  +  A,-(  17  (£+))] 
is  nonzero  for  i  ^  k  but  vanishes  for  i  =  k. 

Next,  using  (8.29  -  8.31)  and  the  Taylor  expansion  of  A*,  we  see  that 
A*(ff(£  +  h-))  -  Afc(Z7(£+))  =  h 

(8-44) 

=  VAjt(£7(£+))  •  (t/(£  +  h-)  -  d(£+))  +  0(|u;(h)|2) . 

If  we  set  jk  =  VAfc(f7(£+))  •  rk(U(£+)),  jk  ^  0  by  hypothesis,  and  use  (8.44),  (8.41)  and  relations 
(8.43)  for  i  ^  k,  we  arrive  at  the  estimate 


jfcWfc(h)  -  h  = 


(8.45) 


O  (Kh)|2)  +  O ( £  |w(s)|  ds)+0{h  |u>(h)|) . 


Adding  (8.43)  for  i  ^  k  with  (8.45)  gives 


■  =  I jkUk(h)  -  h\  +  E  |w;(/i)| 

itk 


(8.46) 


=  O  ( (Hh)|  +  h)  \u(h)\  )  +  0{£  Ms)|  ds ) 

=  0  ( (|«(fc)l  +  h )  <p(h)  )+0{£  <p(s)  ds )  +  0(h2 ) . 
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Since  u >(h)  — ►  0  as  h  — ►  0+,  we  can  choose  £  sufficiently  small  so  that  for  0  <  h  <  6 


(8.47) 


tp(h)  <Ch2+C 


ds 


The  integral  inequality,  in  turn,  yields 

(8.48)  0  <  cp(h)  <C'h 2  for  0  <  h  <  6 


and  thus 


(8.49) 


lim  u>t(h)  _  q  for  i  ^  k , 
fc- o+  h 


lim 

fc-*0+ 


h  Jjt  ' 


This  shows  (8.38).  The  proof  of  part  (ii)  is  virtually  identical,  g 

Proposition  8.5  implies  that  U  has  right  and  left  derivatives  at  any  point  £  which  is  not  an 
accumulation  point  of  <5.  If  such  a  point  £  belongs  to  C  then  U  is  Lipshitz  there,  and  if,  in  addition, 
it  is  an  interior  point  of  I\k  then  /  is  differentiable  there.  It  also  completes  the  picture  regarding 
the  structure  of  the  wave  fans.  We  distinguish  the  following  cases: 

(i)  If  I\k  consists  of  a  single  point  then  the  solution  is  a  shock  wave  satisfying  the  weak  form  of 
the  Lax  shock  conditions  (8.30). 

(ii)  If  I\k  is  a  full  interval  of  points  in  C  the  solution  is  a  k-rarefaction  wave  (provided  that 

VA*  •  r*;  0  on  I\k  which  is  anyway  necessary  for  rarefactions). 

(iii)  In  general  I\k  consists  of  an  alternating  sequence  of  shock  waves  and  k-rarefaction  waves  such 
that  each  shock  adjacent  to  a  rarefaction  from  one  side  is  a  contact  discontinuity  on  that  side. 


* 
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9.  Self-similar  viscous  limits  and  shock  profiles 

In  this  section  we  discuss  the  relation  between  self-similar  viscous  limits  and  shock  profiles  for 
strictly  hyperbolic  systems.  It  was  conjectured  by  Dafermos  [D2]  and  Tupciev  [TU2],  and  proved 
for  systems  of  two  equations  [D2]  that  self-similar  viscous  limits  have  the  internal  structure  of 
traveling  wave  solutions.  We  pursue  here  the  question  in  the  context  of  general  systems. 

Let  £  be  a  point  of  discontinuity  of  U  and  note  that  U(£i:)  satisfy  the  Rankine-Hugoniot 
conditions  (8.3).  Consider  a  sequence  of  points  with  the  property  f  as  e  — ►  0.  Define 

the  function 

(9.1)  Ve(0  =  U£(te  +  eQ,  -00  <  C  <  00  . 

This  introduces  a  stretching  of  the  independent  variable  centered  around  the  point  a  shift  of 
the  shock  speed  £.  The  uniform  estimates  (C*,),  ( S )  imply  that  Ve  is  uniformly  bounded  and  that 

(9.2)  TVC  Fc(-)  =  TVq  t4(6  +e-)  =  TVi  Ue(-)  <  C  . 

Using  Helly’s  theorem  and  a  diagonal  argument  we  establish  the  existence  of  a  subsequence  and  a 
function  V  such  that 

(9.3)  Ue((€  +  e  C)  ->  V(C)  pointwise  for  -  00  <  (  <  00  . 

Proposition  9.1.  Let  £  6  S  and  suppose  that  {£e}  is  a  sequence  of  points  with  — ►  £.  Then  the 
function  U(£)  defined  in  (9.3)  is  continuously  differentiable  and  satisfies  on  (—00,00)  the  traveling 
wave  equations 

dV 

(9.4)  -( [  V  -  V((~) )  +  [  F(V)  -  F (P(f-))  ]  = 
with  initial  conditions 

(9.5)  F(0)  =  lim  tfe(&)  . 

The  limits  lim  F(£)  =:  V±  exist,  are  finite,  and  V+,  V _  satisfy  the  algebraic  equations 

C— +±00 

(9.6)  -t[V-  U(t~) ]  +  [ F(V)  -  F(U(£—)) }  =  0 . 


Proof.  We  evaluate  (8.8)  between  the  points  ££+£C  and  0  and  then  integrate  the  resulting  equation 
in  6  between  £  and  f  +  6,  for  some  6  0,  to  arrive  at 

1  f^+s 

[  -(&  +  £  QU'it,  +  €()+  HUeite  +  £  0)  ]  "  J  J  [  -9U.(9)  +  F(U.(0))  ]  d6 

1  /•{+«  r^e+cC  J  J  /*f+5 

+  jJ  Je  U'(t)  dr d6  =  ±  {Ufa  +  eQ)  -  e-  jf  U&)  d9  . 
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After  an  integration  in  £  we  get 


(9.8) 


r  C  1  [t+s 

Jq  [  - (6 + + «) + F(U'(te + «)) ]  *  -  c  *  y  [~dUe{e) + F^9))] de 

+  -g  Jo  J  Jg  Ue{r)  drdOds  =  tf.(f.  +  £()  -  Ue((e)  --fj'  U'£(9)  dd  . 


j  yC  /•£+*  1‘it+es 

k 

Letting  e  — ►  0  and  using  (9.3),  ( Cb ),  (8.1)  and  (8.10)  we  deduce 


(9.9) 


J\  ~t  V(s)  +  F(V(s ))  ]  ds  -  C  j  ^+'[  -9U(9)  +  F(£f(0))  ]  <M 


1  rZ+S  rZ 

+  C^  ^  Je  U(r)  drd9  =  V(Q  —  V (0) 

From  (9.9),  by  letting  consecutively  6  — *  0+  and  6  — >■  0-,  we  obtain 

(9.10)  J‘  [  -  -  f/(f±))  +  F(VM)  -  F(£l«±))]  ds  =  no  -  no) . 

It  follows  from  (9.10)  that  V(£)  is  a  continuously  differentiable  function  that  satisfies  the 
traveling  wave  equations  (9.4)  and  the  initial  conditions  (9.5).  Since  V  is  of  bounded  variation  on 
IR,  the  limits  lim  V(£)  =:  V±  exist  and  are  finite.  Also,  for  any  integer  n 

C— *±oo 

(9.11)  £+1  [  -  f  (V(s)  -  U((-))  +  F(V(s ))  -  F(U(t-))  ]  ds  =  V(n  +  1)  -  V(n) 

Taking  the  i-th  component  of  (9.11)  and  using  the  mean  value  theorem,  we  deduce  there  are  t\ 
with  n  <  txn  <  n  +  1  such  that 


(9.12)  -tiV'ti)  -  U*(t-))  +  F^VifJ)  -  Fi(U((-))  =  V\n  +  1)  -  V‘(n) ,  i  =  1 ,  ... ,  N 


Letting  n  — ♦  oo  shows  that  V+  is  an  equilibrium  for  (9.4).  Similarly,  V_  satisfies  (9.6).  | 

The  function  V  as  well  as  the  limiting  values  V±  depends  on  the  choice  of  the  sequence  {£e}. 
For  several  choices  of  {£e}  it  may  happen  that  the  traveling  wave  disintegrates  to  a  constant 
solution.  Two  questions  arise:  (i)  Is  it  always  possible  to  choose  {fe}  so  that  the  resulting  V  does 
not  disintegrate  to  a  constant  solution  of  (9.4).  (ii)  What  is  the  relation  of  (7(£— ),  U (£+)  and 
nontrivial  heteroclinic  orbits. 

Proposition  9.2.  Let  £  6  S  be  fixed  and  suppose  the  set  of  solutions  to  (9.6)  is  not  connected. 
There  exists  a  sequence  of  shock  shifts  {£e}  such  that  the  resulting  V  in  (9.3)  is  a  nontrivial 
heteroclinic  (or  homoclinic )  orbit. 
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Proof.  Suppose  the  solution  set  of  (9.6)  with  £  fixed  is  contained  in  two  open  sets  O-  3  {/(£—) 
and  Q+  9  Z7(£+)  with  O _  D  0+  =  0.  Because  of  (Cb)  we  may  restrict  attention  to  a  ball  Bm 
containing  Ue.  For  a  large  integer  n  it  is  {/(£  —  ■-)  €  O _  and  £f(£  +  ~)  £  Choose  £n  such 
that  £4B(f  -  £)  €  £>_,  £1£b(£  +  £)  €  0+.  There  exists  {&„}  satisfying  £-£<&„<£+£ 
and  U£n({£n)  €  BM  ~  ( O _  U  <9+).  Along  a  subsequence  (£n  -»■  £  and  U£n((£n)  ->  ^(0)  with 
F(0)  ^  O-  U  O+ .  The  resulting  V  is  a  nonconstant  solution  of  (9.4)  connecting  two  equilibria  V_ 
and  V+.  | 

The  hypothesis  in  Proposition  9.2  is  violated  only  for  shocks  associated  with  a  linearly  degen¬ 
erate  characteristic  field  :  VA^(17)  •  r*(£l)  =  0  for  all  U  (c./.  Section  10)  .  Addressing  (ii)  is  quite 
complicated  at  the  full  level  of  generality.  We  give  one  result  indicating  what  can  happen  if  there 
is  a  finite  number  of  equilibria  in  Bm,  the  range  of  variation  of  U£. 

Proposition  9.3.  Let  £  £  <S  and  suppose  that  (9.6)  has  a  finite  number  of  solutions  in  Bm-  There 
exists  a  subsequence  en  — *  0  and  choices  {fi£n},  {^2e„}  °f  the  shock  shifts  such  that  £i£b  <  ^2ea, 
£le„  -*•  £>  be„  -*■  Z, 

(9.14)  UCu(t Un  +  SnC)  -  MO ,  V£%(t2.n  +  £nO  -  MO ,  pointwise  for  -  oo  <  C  <  oo  , 

and  the  resulting  V\,  V-i  are  solutions  of  (9.4)  that  satisfy  M— oo)  =  £!(£—),  M+00)  =  U(£+)- 

Proof.  Let  Bm  be  the  ball  where  the  solutions  U£  range,  and  suppose  (9.6)  has  a  finite  number 
of  solutions  £/(£—),  £/(£+)  and  U\, ... ,  Uj.  Fix  two  open  balls  B-,  B+  and  an  open  set  O  with 
the  properties:  B- ,  B+  and  O  lie  inside  Bm,  B-  is  centered  at  U (£-),  B+  is  centered  at  £/(£+), 
O  contains  U\,  ... ,  Uj,  and  the  distances  among  any  two  of  the  sets  B-,  B+  and  O  are  strictly 
positive.  Since  U  is  of  bounded  variation,  we  can  fix  S  >  0  such  that  U(6)  €  B _  for  0  £  [£  -  6,  £) 
and  tf(0)  €  5+ for  0  £(£,£  +  £]. 

Consider  a  convergent  sequence  :  UCn  -*  U  pointwise  on  1R.  In  the  sequel  we  will  be  extracting 
appropriate  subsequences  that  are  denoted  again  by  U£n.  Choose  no  such  that  UCn(£  -  6)  £  B-, 
UCn(£  +  0)  £  B+  for  n  >  no-  For  each  n  >  no  choose  points  a 'n,  Aln ,  b'n,  B'n,  i  =  1,  ... ,  K(n),  in 
the  interval  Is  =  [£  —  0,  £  -f  0]  to  be  respectively:  a\  =  £  —  0,  6*  the  first  point  where  U£n  enters 
the  ball  i?+,  A\  the  last  point  in  (a*,60  that  U£n  exits  2?_,  a2n  the  first  point  after  bln  that  U£n 
enters  the  ball  J5_  (if  applicable),  B\  the  last  point  after  b\  that  U£n  exits  B+,  and  so  on  until 
finally  Bn  ^  =  £  +  0.  These  are  defined  as  follows: 

&n  =  inf {0  >  a'n  :  U£n(0)  £  B+ }  ,  A'n  =  sup{0  £  (a'n,b'n)  ■  U£ .(0)  €  B.}  , 
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°n+1  =  inf{0  >  K  :  U£n(0)  eB-},  B'n  =  sup{0  €  (^,<+1) :  UCn(0)  €  B+}  ; 

if  in  the  i-th  step  ajj1"1  is  not  well  defined  then  i  =  K(n )  and  B'n  =  £  +  6.  Since  UCn  is  of  uniformly 
bounded  variation,  it  can  go  back  and  forth  between  to  B+  at  most  a  finite  number  of  times, 
thus  K{n )  is  bounded.  By  restricting  to  subsequences  we  may  assume  K{n)  is  some  positive  integer 
K  for  large  n,  and  that  a'n  — ►  a\  A'n  — ►  A*,  b'n  — ►  b *  and  B'n  -*  B',  t  =  1,  ... ,  K,  as  n  -*■  oo. 

By  construction  <  <  K  <  K  <  B'n  <  a{j+1  and  UCn  has  the  behavior  UCn  (6)  €  Bm  -  B+  on 
(a'^A'n),  USn(0)  €Bm-  ( B _  U  B+)  on  (A'n,b'n)  and  UCn(0)  €  BM  -  B-  on  ( b'n,B'n ).  As  a  result 
the  limits  a*,  A *,  6'  B'  have  the  following  properties:  (i)  a'  <  A*  =  b'  <  B'  =  a,+1,  (ii)  if  Bl  <  £ 
then  A'  =  6*  =  Bl  =  a,+1,  and  (iii)  if  a'  >  £  then  B'~l  =  a*  =  A*  =  b\  To  see  (ii)  suppose  that 
Bx  <  if  A'  <  a*+1  there  is  0  <  £  such  that  UCn{0 )  £  for  large  n.  Passing  to  the  limit  we  see 
that  U{0)  &  B-  a  contradiction.  The  rest  are  proved  by  similar  arguments. 

In  what  follows  we  fix  j  to  be  the  first  index  that  B3  =  £  and  k  to  be  the  last  index  that 
ak  =  £.  Then  we  have  the  ordering 

(9.15)  a3  <  A3  =  b3  =  B3  =  ...  =  |  =  ...  =  ak  =  Ak  =  bk  <Bk 
for  any  index  between  j  and  k . 

*  Consider  first  the  case  that  (9.6)  has  precisely  two  solutions  {/(£—)  and  U(£+).  Set 

fie.  =4;,  Flen(C)=£U^n+£nC), 

(9.16) 

6..  =bk,  V2Cn(0  =  UCn{bk  +  en0. 

Then  VieB(0)  lies  on  dB-  and  V2C„(0)  lies  on  dB+.  Along  a  subsequence  V\Cn  and  V2Cn  converge 
pointwise  to  a  solution  of  (9.4)  and  the  limits  V,(±oo)  =  V{±  exist  and  are  finite.  Since  no  solutions 
of  (9.6)  lie  on  the  boundaries  of  j B_  and  B+  the  resulting  traveling  waves  are  nontrivial.  From  the 
definition  of  V\En  and  V2Cn  follows  that 

(9.17)  Vu9(0  ?B+  for  <—*k  <C<0,  V2e„(C)  i  B.  for  0  <  C  <  ^  . 

£n  £n 

Since  lim  a3n  =  a3  <  £  =  lim  A3n  and  lim  b„  =  £  <  Bk  =  lim  Bk, 

(9.18)  Vi(C)  $  B+  for  -  oo  <  C  <  0 ,  V2(C)  t  B-  for  0  <  C  <  oo  . 

As  {/(£-)  is  the  only  equilibrium  in  Bm  -  B+  and  U(£+)  the  only  equilibrium  in  Bm  -  -B-  it 
follows  Vj(-oo)  =  17(f-)  and  F2(+oo)  =  f/(£+). 

Suppose  next  that  (9.4)  has  more  than  two  equilibria.  If  U£n  never  enters  O  the  previous  proof 
shows  the  desired  result.  If  UCn  enters  O,  we  restrict  attention  to  the  interval  [a^,  b3n]  and  proceed 
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as  Mows:  Note  that  U£n(0)  <E  BM  -  B+  on  «, b>n)  and  U£„(b}n)  €  dB+.  As  in  the  previous  step 
we  choose  points  c'n  <  Cln  <  d'n  <  D'n,  i  =  1, ... ,  K(n),  with  the  properties:  c*  =  aJn,  d\  the  first 
point  after  c\  that  U£n  enters  O,  C*  the  last  point  before  d\  that  U£n  exits  £?_;  if  U£n  reenters  5_ 
then  we  define  c\  to  be  the  first  point  after  d\  that  U£n  enters  the  ball  B-,  define  D\  to  be  the  last 
point  before  c\  that  U£n  exits  O,  and  reiterate  the  above  procedure;  if  U£n  does  not  reenter  B- 

then  we  set  D\  to  be  the  last  point  of  exit  from  O  before  touching  dB+  and  stop  at  this  step.  As 

the  sequence  { U£n  }  is  of  uniformly  bounded  variation  the  process  will  conclude  in  a  finite  number 
of  steps.  By  restricting  to  subsequences  we  may  assume  K(n)  =  K  <  oo,  cxn  — *■  c*,  C'n  — *•  C', 
d'n  -»■  d',  D'n  -*  D'.  Again  if  d'  <  {  for  some  i  then  C'  =  d'  =  D'  =  c,+1.  Let  l  be  the  first  index 
such  that  Dl  =  e.  Then  Dl~l  =  c‘  <  Cl  =  d‘  =  Dl  =  £.  If  we  set 

(9.W)  (u.  =  C‘n ,  V„. (0  =  V..  K  +  £nC) 

then  V\Cn  satisfies 

(9.20)  VUn  (C)  i  O  U  B+  for  C"  -  <  C  <  0 

and  the  resulting  traveling  wave  V\  hais  the  property  Vi(— oo)  =  A  similar  construction 

shows  the  second  part  of  the  Proposition.  | 

Proposition  9.3  shows  that  if  £  is  a  point  of  discontinuity  of  a  solution  U  arising  via  self-similar 
viscous  limits  then  there  exists  one  heteroclinic  orbit  of  (9.4)  that  emanates  from  {/(£—)  and  one 
that  concludes  to  {/(£+).  ^  expected  that  in  general  this  will  be  the  same  heteroclinic  orbit. 
However,  if  more  than  two  states  in  Bm  satisfy  the  Rankine-Hugoniot  conditions  (9.6)  at  a  given 
£  £  iS,  or  if  multiple  heteroclinic  connections  between  two  equilibria  are  possible,  then  the  precise 
relation  between  self-similar  limits  and  shock  profiles  requires  a  detailed  analysis  of  the  heteroclinic 
orbits  (the  proof  is  suggestive  as  to  what  possibilities  must  be  excluded).  In  specific  examples  it 
usually  happens  that  there  is  a  single  shock  profile  connecting  U(£— )  to  {/(£+).  ^  however 
possible  that  there  are  intermediate  states  Vjy  j  =  1,...,  J,  satisfying  (9.6)  and  a  chain  of  shock 
profiles  with  the  same  shock  speed  £  that  connect  succesively  t7(f—)  to  Vj,  each  of  the  points 
Vj  to  the  next,  and  Vj  to  U(£+).  The  latter  situation  occurs  for  the  equations  of  elasticity  in 
the  presence  of  multiple  inflection  points  in  the  stress-strain  relation,  for  specific  positions  of  the 
Riemann  data  relative  to  the  stress-strain  curve  [Tz2]. 
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10.  Comparisons  with  the  classical  solution  of  the  Riemann  problem 

In  this  section  we  compare  the  classical  solution  of  the  Riemann  problem  with  the  solution 
obtained  via  self-similar  viscous  limits.  For  systems  of  strictly  hyperbolic  conservation  laws  the 
classical  solution  of  the  Riemann  problem  is  based  on  a  detailed  study  of  elementary  solutions  of 
rarefaction  waves  and  shock  waves,  and  was  established,  for  \U+  —  U-\  small,  by  Lax  [Lai]  in  the 
genuinely  nonlinear  case  and  by  Liu  [Lii,  L4]  in  the  general  case. 

Fix  Uq.  Let  TZk  =  TZ^Uq)  be  the  integral  curves  of  the  vector  field  r*  emanating  from  Uq. 
Rarefaction  wave  solutions  take  values  on  the  curves  TZk-  Shock  waves  emerge  by  solving  the 
Rankine  Hugoniot  conditions 

(10.1)  s{U-U0)  =  F(U)-F(U0). 

For  U  near  Uq ,  the  set  of  solutions  of  (10.1)  consists  of  N  smooth  curves  Sk  =  Sk(Uo)  tangent  to 
TZk(Uo)  at  Uq ,  k  =  1,  ... ,  N.  Each  Sk  is  associated  with  the  fc-th  characteristic  field,  it  is  defined 
by  parametric  equations  U  =  £4(r)  and  s  =  Sk(r)  for  |r|  small,  and  the  parametrization  may  be 
»  arranged  so  that 

£4(0)  =  U0 ,  U(0)  =  rk(U0), 

*  i 

(10.2)  sk(0)  =  Xk(Uo) ,  ifc(0)  =  -  VA*(0b)  •  rk(U0) , 

Xk-i(Uk(T))<»k(r)<Xk+1(Uk(T)). 

A  state  Uk(r)  €  Sk(Uo)  gives  rise  to  a  shock  wave  solution  with  speed  Sfc(r),  left  state  Uq,  and  right 
state  Ukij).  Liu  [L4]  performed  a  detailed  study  of  the  shock  curves  and  proposed  the  following 
shock  admissibility  criterion.  A  shock  {Uq,  Uk{r),  ■Sfc(r))  is  admissible  if  it  satisfies 

(E)  »*(r)  <  Sk{t )  for  t  between  0  and  r . 

Using  ( E )  Liu  obtained  a  unique  solution  the  Riemann  problem. 

Consider  the  solution  U  constructed  via  self-similar  viscous  limits  in  the  previous  sections. 
U(£)  takes  values  in  a  small  ball  B^{U-),  the  wave  speeds  X k{U{£))  are  separated,  and  U{£)  has 
the  properties  indicated  at  the  end  of  Section  8.  Each  wave  fan  is  studied  separately  and  we 
distinguish  three  cases: 

(i)  At  is  genuinely  nonlinear :  VA k(U)  •  rk(U )  ^  0  for  all  U . 

«  For  a  genuinely  nonlinear  characteristic  field  the  shock  speed  Sk(r)  is  increasing  in  one  direction 

of  the  shock  curve  Sk(Uo)  and  decreasing  in  the  opposite  direction.  Contact  discontinuities  are 
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excluded  for  weak  shocks.  The  behavior  of  U  on  I\k  simplifies  considerably:  Either  I\k  is  empty, 
or  I\k  consists  of  a  single  point  of  jump  discontinuity  £  and  U  satisfies  the  Lax  shock  conditions 

(10.3)  A*(£f(£+))<£<A*(lf(£-)), 

or  I\k  is  a  full  interval  of  points  of  continuity  and  the  solution  is  a  fc-rarefaction  wave  on  I\k . 
Therefore,  for  genuinely  nonlinear  and  strictly  hyperbolic  systems,  the  emerging  structure  of  U  is 
identical  to  that  determined  by  Lax  [Lai]. 

(ii)  A k  is  linearly  degenerate  :  VXk(U)  •  Tk{U)  =  0  for  all  U. 

For  a  linearly  degenerate  characteristic  field  the  k- th  shock  curve  emanating  from  Uq  is  given  by 
U  =  s  =  sk(r )  where 

(10.4)  Sfc(r)  =  Xk(U0)  ,  ~j~~(r)  =  rk(Uk(r))  ,  Uk( 0)  =  U0 . 

A  version  of  the  converse  is  also  true:  If  (10.1)  has  a  curve  of  solutions  U{r)  corresponding  to 
s(r)  =  Sq  fixed  then  U(r)  =  rfc(/7(r)),  s0  =  A k(U(r ))  for  some  fc,  and  the  k- th  field  is  linearly 
degenerate.  Since  A*  remains  constant  on  the  curves  rarefaction  wave  solutions  are  not  possible 
for  linearly  degenerate  characteristic  fields.  A  close  look  in  the  proof  of  Propositions  8.4  and  8.5 
shows  that  it  is  not  possible  that  hk  is  a  full  interval.  Therefore,  either  I\k  is  empty,  or  it  consists 
of  a  single  point  of  jump  discontinuity  and  U  is  a  contact  discontinuity. 

(iii)  The  curves  TZk  intersect  the  set  { U  :  VXk(U)  -  rk(U )  =  0}  at  discrete  points. 

The  solution  U  cannot  be  further  simplified  in  this  case.  The  relation  with  the  Liu  shock  ad¬ 
missibility  criterion  ( E )  is  established  indirectly,  using  Proposition  9.3  on  the  relation  between 
self-similar  limits  and  shock  profiles,  in  conjunction  with  results  of  Liu  [L^],  Majda  and  Pego  [MP] 
on  the  relation  between  shock  profiles  and  (a  strict  inequality  version  of)  condition  ( E ).  Majda 
and  Pego  [MP,  Theorem  3.1]  prove  that  given  two  states  !/(£—)  and  17(£+)  in  a  small  ball 
satisfying  the  Rankine-Hugoniot  conditions  for  some  speed  £,  a  shock  profile  lying  in  B^(U-)  ex¬ 
ists  connecting  U (£— )  and  U (£+)  if  and  only  if  condition  (E)  is  satisfied  as  a  strict  inequality. 
Moreover,  there  exists  at  most  one  trajectory  F(£)  of  (9.4)  connecting  U (£-)  and  U (£+)  which 
remains  in  B^U-)  for  all  (. 

Fix  £  6  S  n  I\k  and  consider  the  set  of  all  solutions  to  the  Rankine-Hugoniot  conditions  that 
are  compatible  with  (8.30).  If  U{i~)  and  U (£+)  are  the  only  states  with  this  property  then  there 
is  a  shock  profile  connecting  them  and  the  shock  speed  £  satisfies  the  strict  condition  ( E ).  If  there 
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are  more  than  two  such  solutions  of  (9.6)  then  there  is  a  shock  profile  in  B^U-)  connecting  U (f — ) 
to  some  state  V:  and  another  shock  profile  connecting  a  state  V,  to  {/(£+).  It  is  expected  that  in 
this  case  there  will  be  a  chain  of  shock  profiles  that  connect  {/(£—)  through  intermediate  states 
with  (eventually)  U(£+). 
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